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Abstract. We offer connections between some problems of PDE, geometry, 
algebra, analysis and physics. The uniqueness of the solution of the Dirichlet 
problem and some another boundary value problems for the string equation 
inside of an arbitrary biquadratic algebraic curve is considered. It is shown 
that the solution is non-unique if and only if corresponding the Poncclet prob- 
lem for two conies has a periodic trajectory. Similarly some other problems are 
proved to be equivalent to it. Among them there are the solvability problem of 
the algebraic Pell-Abel equation and an indeterminacy problem of a moment 
problem that generalizes well-known trigonometrical moment problem. Solv- 
ability criterions of above problems can be presented in form 9 € Q where the 
number 8 = m/n is connected with the concrete problem data. We demon- 
strate also an intimate relation of the above mentioned problems with such 
problems of the modern mathematical physics as elliptic solutions of the Toda 
chain, static solutions of the classical Heisenberg XY-chain and biorthogonal 
rational functions on elliptic grids in the theory of the Pade interpolation. 
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1. Introduction 

This work is devoted to new connections between some problems of mathematics 
such as some ill-posed boundary value problems in a bounded semialgebraic domain 
for partial differential equations, a moment problem, the Poncelet problem of the 
projective geometry, the algebraic Pell-Abel equation and some other problems, 
recently revealed by authors. 

Study of ill-posed boundary value problems in bounded domains for partial differ- 
ential equations go back to J.Hadamard and they are a popular object of modern 
infcstigations (see s. I2.1j) . In this paper we will examine generally the Dirichlet 
problem for the string equation 

u xy = in fi, (1.1) 

u\ c = 4> on C = dCl, (1.2) 

the solution uniqueness for which is the problem of existence of a nontrivial solution 
of the homogeneous Dirichlet problem 

u\ c = on C. (1.3) 

The functions u, <fi are assumed to be complex- valued functions of real variables. 
We will consider this and other boundary value problems in semialgebraic domains, 
the boundary of which is given by so-called biquadratic algebraic curve 

2 

F(x,y) := a lk x l y k = (1.4) 

i,fe=0 

where x i , y k are powers. We will consider canonical forms of the curve (|1.4| to 
which the generic curve can be transformed by linear-fractional replacements and 
we will give criterions of uniqueness breakdown in the form 

reQ (1.5) 

where the number r is determined by the curve C . Our investigations are based 
on a study of the John mapping generated on C by characteristics of the equation 
(jl.lj) . see propositions[29](s. [6]) and[T9l(s. 13. 5p where we use that this John mapping 
becomes an usual shift after a transform on universal covering group of the variety 

C3J). 
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We will have observed a remarkable connection of this problem with the Poncelet 
problem. The Poncelet problem is one of famous problems of projective geometry 
and it by itself has numerous links with a set of different problems of analysis 
and physics (see [13], [10] [26] and below s. 14. ip . We will prove that for generic 
biquadratic curve C the Dirichlet problem has a non-unique solution if and only 
if corresponding the Poncelet problem has a periodic trajectory (proposition 1251 s. 
14. 5p . And therefore we will give a new criterion of periodicity in the form (|1.5|) that 
differs from the well-known Cayley criterion (see s. 14. 6[) . 

Remind that the existence of a periodic trajectory in the Poncelet problem means 
that each trajectory is periodic by the big Poncelet theorem [15] . Note that different 
cases of disposition of conies give different cases of curves C and a new setting of 
the Dirichlet problem (11.1|) . ()1.3|) for unbounded domains is more suitable than the 
classical setting (see s. I2.2|) . 

We will observe one more remarkable connection of the Poncelet problem with 
the solvability problem of algerbraic the Pell- Abel equation 

P{tf - R{t)Q{tf = L (1.6) 

where for given polynomial R(t) of one variable t one should find polynomials P, Q 
and a constant L satisfying the equation. The Pell-Abel equation solvability is 
one of famous problems also and this algebraic problem by itself has numerous 
connections with many problems of analysis (see below s. [5}. 

We will prove (see proposition l27l of s. 15. 2p that the Pell- Abel equation (|5.ip 
with a polynomial R of the third or fourth order has a solution if and only if 
corresponding the Poncelet problem has a periodic trajectory with an even period. 
On this way we obtain a new criterion in the form (|1.5p which differs from well- 
known others, for example, from well-known Zolotarev's porcupine (see [43]). 

We will show that the same condition (|1.5p is a criterion of solution uniqueness 
to within an additive constant of the Neumann problem 

u*.\c = ip (1-7) 

for the same equation (II. ip in the same domain where w„, is a derivative with 
respect to the conormal (statement [3] of s. 12. 3p . We will show also that the 
same condition is a criterion of indeterminacy of a moment problem on the curve 
C (statement [5] of s. 12.31) . By means of the duality equation-domain we will 
obtain an equivalent problem (|2.22p in the form of a hyperbolic (in some cases 
of corresponding Poncelet problem) equation of the fourth order with only two 
boundary data on characteristics instead four as in the boundary value problem of 
the Goursat type, where we will have solution uniqueness almost always. 

We will have observed some connections with a problem on classical XY-spin 
chains and a problem from the theory of Toda chain. We will have shown an equiv- 
alence of considered problems and also we give an interpretation of this criterion in 
terms of the John mapping. About links with some other problems of analysis see 

[SH], in], Ha], m]. 

Note that some results of present work have published already as short and 
incomplete fragments in papers [IS], [IS], [213] , 

Remark also that some explicit necessary and sufficient conditions of uniqueness 
breakdown of Dirichlet problem (and some others boundary value problems) solu- 
tion for partial differential equations with constant coefficients have obtained (see 
e.g., [16] and [17]) in an arbitrary ellipse. Answers in that works were formulated 
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just in the form of condition (|1.5j) and prompted answers in our present investiga- 
tions. In this paper we would like to open a way to examine ill-posed boundary 
value problems for partial differential equations in some more complicated domains 
than a circle along with a survey of another fields of mathematics about topics with 
equivalent contents. 

The main goal of the present paper is thus to show numerous relations of the con- 
sidered problem with remote branches of mathematics and mathematical physics. 
The paper is organized as follows. 

In Section 2 we describe a theory of uniqueness of Dirichlet problem and the 
John algorithm. In section 3 we apply this technique to the concrete choice of the 
biquadratic curve. In Section 4 we consider the solution of the Poncelet problem and 
show its relation with the John algorithm for the generic biquadratic curve. Section 
5 is devoted to the Abel problem of reducing of elliptic integrals to elementary 
ones. We show how this problem can be formulated in terms of the Poncelet 
problem (or the John algorithm for the biquadratic curve). In Section 6 we consider 
relation of our problems with Ritt's problem of existence of periodic functions with 
a nontrivial multiplication property. Finally, in Section 7 we consider 3 problems in 
mathematical physics which are related with the Poncelet (or the John algorithm) 
problem: static solutions of the classical XY Heisenberg chain, elliptic solutions of 
the Toda chain and elliptic grids for biorthogonal rational functions in the theory 
of rational Pade interpolation . 



2. Boundary value problems in a domain for the string equation 

2.1. Bibliographical remarks. Investigations of ill-posed boundary value prob- 
lems in bounded domains for partial differential equations go back to J.Hadamard 
[30] and then A.Huber [32] who for the first time noted nonuniqueness of the solu- 
tion of the Dirichlet problem for the equation of string vibration (string equation) 
in a rectangular. Boundary value problems in bounded domains for nonelliptic par- 
tial differential equations were regularly studied on the whole in a parallelepiped, 
besides in domains with a general boundary questions of solution uniqueness of the 
Dirichlet problem for the hyperbolic equation of the second order on the plane (see 
reviews and results in [48] were usually studied. In work [14] D.Burgin and R.Duffin 
have examined the homogeneous Dirichlet problem for the equation u t t — u xx — 
in a rectangular {0<t<T;0<x< X}. It is shown that if the ratio T/X is irra- 
tional, uniqueness in space of continuously differentiated functions with summable 
second derivatives takes place. Theorems of existence of solutions in classical spaces 
are established, and smoothness of the solution is that more than smoothness of 
boundary function is big and than the number T/X is worse approximated by ratio- 
nal numbers. The Neumann problem is considered also. In works by B.Yo.Ptashnik 
and his pupils boundary value problems in a parallelepiped for a wide class of the dif- 
ferential equations and systems are investigated, see [48]. All these works (excepting 
the mentioned work [6]) are based on the methods, essentially using representation 
of domain as a topological product. 

For nonrectangular domains the Dirichlet problem for the string equation was 
studied in connection with a number of the Denjoy-Poincare rotation (see, e.g., 
Z.Nitetsky's book [46] ) of a homeomorphism of domain boundary, constructed on 
characteristics of equation (so-called an automorphism of characteristic billiards by 
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Fritz John [34]). A connection of properties of the Dirichlet problem with prop- 
erties of this homeomorphism has used even in mentioned works by J.Hadamard 
and A.Huber. To the analysis this connection has undergone in F.John's work [34j . 
In works by R.A.Alexandrjan and his pupils investigations of this problem and, 
in particular, of this connection was continued ((3], [4], [47], [5]). The question on 
uniqueness of the solution of the Dirichlet problem in this ideology for domains 
that are convex with respect to characteristics families, should be transformed to a 
question on irrationality of number of rotation or, that is the same, to a question on 
presence a continuum set of finite orbits (cycles) of the discrete dynamic system gen- 
erated by mentioned John homeomorphism. The same questions in connection with 
an asymptotic behavior of the solution of the Sobolev equation that describe surface 
oscillations of a fluid filling a body which is flying in atmosphere, were investigated 
by the Siberian mathematicians T.I.Zelenjak, I.V.Fokin and others ([62], [21]). Re- 
searches of the string equation are included in well-known Yu.M.Berezansky's book 
[12] also, they give possibilities to build domains with angles in which the homoge- 
neous Dirichlet problem is weakly solvable and well-posed concerning the right part 
and small moves of boundary of the domain, leaving angles in specified limits. Note 
more that the case of an ellipse was considered in works by A. Huber [32] . Alexan- 
drjan [3], V.I.Arnold [5], for small smoothness and for more general equations see 
the book [T7j. 

2.2. John condition. For the problem (jl.ip . (|1.2[) in some general bounded do- 
main Fritz John [34] considered a remarkable transformation T : C — > C of 
Jordan boundary into itself, allowing to do some conclusions about properties of 
the Dirichlet problem in Q. Let us describe it. 

Let n be arbitrary bounded domain which is convex with respect to character- 
istics of the equation i.e. it has the boundary C intersected in at most two 
points by each straight line that is parallel to x- or y-axes. We start from arbitrary 
point Mi on C and consider a vertical line passing through M\. Generally, there are 
two points of intersection with the curve C: Mi and some point M 2 . We denote I\ 
an involution which transform M\ into ■ Then, starting from , we consider a 
horizontal line passing through M^. Let M3 be the second point of intersection with 
the curve C. Let I2 be corresponding involution: I2M2 = M3. We then repeat this 
process, applying step-by-step involutions I\ and I2. Denote T = L 2 Ii,T~ 1 = 
This transformation T : C — ► C gives us a discrete dynamical system on C, i.e. an 
action of group Z and each point M £ C generates an orbit {T n M\n £ Z}. This 
orbit can be finite or denumerable set. The point M with finite orbit is called a 
periodic point and smallest n, for which T n M = M, is called a period of the point 
M. In the paper [34] the uniqueness breakdown in the problem (|1.1[) . (|1.2[) have 
studied in connection with topological properties of the mapping T for the case of 
even mapping T. The mapping T is called to be even or preserving an orientation 
if each positive oriented arc (P, Q) with points P,Q £ C transforms into positive 
oriented arc (TP,TQ). F.John have proved several usefull assertions, among of 
which we extract the following one. 

Sufficient condition of uniqueness. The homogeneous Dirichlet problem ljl.l|> . 
(|1.3|) in the bounded domain has only a trivial solution in the space C 2 (£l) if the 
set of periodic points on C is finite or denumerable. 

It has been selected four possible cases of dynamical system behavior: 

I) all points are periodic (then their periods coincide); 
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II) there are periodic and nonperiodic points; 

III) there are no periodic points and there is not a point, the orbit of which 
{...,T~ 1 P,P,TP,...,T n P,...} is dense in C; 

IV) there are not periodic points and there is a point, the orbit of which 
{...jT-^PiTP, ...,T n P, ...} is dense in C (transitive case). 

In the work [34] it has shown that for a C 2 -smooth curve the case III) is impos- 
sible. For the case II) it has proved that there is an arc Dq on C such that each 
two arcs from D Q , I\D a , TD Q , IiTD , T 2 D Q , ^r^Do, T n D , ... have not any 
common point. Note that for the analytical boundary this is impossible because 
in this case T is a diffeomorphism. In the case IV) it has proved that the Denjoy- 
Poincare rotation £ ([46]) of the John mapping T is irrational and T is topologically 
equivalent to a turn of an unit circle about the angle 7r£ (i.e. there exists a home- 
omorphism h from C onto the unit circle S such that the mapping hTh^ 1 : S — ► S 
is a turn about the angle 7r£). 

For the case when every point of C is periodic it have proved the all periods are 
coincided. But here we don't know anything about solution uniqueness, although 
R.A. Alexandrjan have shown in this case [3] that there is a generalized solution of 
the problem (|l.ip . (|1.2p which can be built as a piecewise constant function. 

We will assume that the domain has boundary C satisfying the condition: 

The curve C is smooth and each characteristic line either doesn't intersect 

the curve C or touchs it at a point, named a vertex, or cuts C at two points. (2.1) 

Note that above John's condition of convexity with respect of characteristic 
directions will be fulfilled under the condition (|2.ip on C = dil in the case of a 
bounded domain f2. 

For the case of bounded domain with a biquadratic boundary we will see that 
above sufficient uniqueness condition is also necessary, moreover, it will be so even 
for cases when the curve C is unbounded but then we should change the setting of 
the problem. Namely, along with the usual setting of the uniqueness property: 

the examined bounded domain is such that the homogeneous Dirichlet 

problem &l.l)) . n~H\) has only trivial solution in the space C 2 (fl) n C(fl), (2.2) 

for cases when the curve C is unbounded we will examine the following modification 
of uniqueness property for the homogeneous Dirichlet problem : 

the examined curve C is such that each analytic in real sense solution 

in R 2 of the equation hl.l}) with the property hl.3\) is only zero solution. (2.3) 

Note, the assumption "analytic" is introduced in order to allow unbounded curves 
C where there exist characteristic lines which do not intersect C and are between 
curve branches. Without this assumption on solutions for such curve and e.g. with 
an assumption of infinite smoothness for solutions one may build a simple example 
of smooth in R 2 nontrivial solution of the equation with the property (|1.3[) . 
We will call a modified setting the problem to find an analytic in real sense solution 
in R 2 of the equation with the property (|1.3[) . 

Now we will give the John's proof of the sufficient condition of uniqueness in 
order to show that John's arguments are valid also for the problem 12.31 

Proposition 1 . Under the condition (|2.1| if the mapping T is transitive acted on C 
then 1) the uniqueness property (|2.2p holds for the usual setting of the homogeneous 
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Dirichlet problem (fTTj) . (fL~3j) in a bounded domain and 2) the uniqueness property 
(|2.3p holds for the modified setting in, possible, unbounded domain. 

Proof. Let a function u € C 2 (Q) n C(f2) be a nontrivial solution of the problem 
(|1.1[) . (|1.3|) in a domain with the condition (|2.ip . As it is well known, there exist 
two functions u\, u 2 of the class C 2 depending on one variable such that u{x, y) = 
U\(x)+u 2 {y) which we will write for any point P £ C as u(P) — u\{P)+u 2 {P). For 
the case of the property (|2.3p consider a domain f2 C M. 2 of points P = (xo, yo) el 2 
for which there exists a pair of different points of intersection C H {x — xq} of the 
curve with corresponding vertical line and also there exists a pair of different points 
of intersection C D {y — yo} of the curve with the horizontal characteristic line. 
Then for any point P G f2 using definitions we easily obtain: u\(I\P) = u\(P), 
u 2 (I 2 P) = u 2 (P); u 2 {P) - u 2 {TP) = u 2 {P) - u 2 (hP) = u(P) - u{hP) = 0; 
Ui(P)-u 1 (T- 1 P) = u l (P)-u l {I 2 P) =u(P)-u(I 2 P) = 0. From what follows that 
equalities U\(P) = Ui(T n P), u 2 (P) — u 2 (T n P) hold for each integer n. Continuity 
gives us that ui(P) — ui(Q), u 2 {P) — u 2 (Q) for any point Q from closure of the 
orbit of P. Because for the transitive acting the closing of the orbit of any point 
coincides with C then u\ = const and u 2 = const and, therefore, u = in O. In 
the case of the property ()2.3|1 the analyticity allows us to continue the zero solution 
onto the plane. □ 

Above we have noted that for any analytic boundary only two cases are admis- 
sible: periodic (I) and transitive (IV by John). Therefore we give the following 
setting of the periodicity problem for John mapping : 

What curve from a given class of curves has the property: 

The John mapping T : C — > C has at least one periodic point. (2-4) 

Then as we have known already, the all points are periodic on our curve. 

Along with above settings we will consider also a setting with complex John 
mapping. If C G C 2 is a biquadratic complex curve (i.e. one-dimensional complex 
variety (|1.4|l ) then it is fulfilled the property of type 12.11 

Almost each "vertical" line x = xq intersects the curve C at two different 

points and almost each "horizontal" line, too. (2-5) 

Let C G C 2 be an analytic curve with the same property. Then we can construct 
a John mapping T in the same way as in real case. And we can ask a similar 
question about the periodicity problem for the complex John mapping : 

What curve from a given class of curves has the property: 

The complex John mapping T : C — > C has at least one periodic point. (2.6) 

Corresponding complex setting of the problem (jl.l|) . l|1.3|) we will use is the fol- 
lowing: 

To find two meromorphic functions f(x), g(y) of one complex variable such 

that f(x) + g(y) = as soon as (x,y) G C C C 2 . (2.7) 

Note that for the case of transitive action of the mapping T the proof of the 
proposition [1] one can easy transfer onto this complex case: 

Proposition 2. For a biquadratic complex curve (|1.4| if the complex John mapping 
T transitive acts on C then the complex problem (|2.7| has only trivial solution. 
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Below we will consider biquadratic curves C satisfying the property (|2.1j) and 
will give an explicit criterion distinguishing the cases of periodicity and transitivity 
and in the first case we will build an explicit nontrivial solution of the problem 
(jl.ip , (II. 3|) in sense of the settings 12.21 or 12.31 and as an intermediate the setting 
(|2.7|) . In the case of transitivity each solution is proved to be zero. 

2.3. Boundary value problems and a moment problem. In this subsection 
we are going to indicate an equivalence of properties of some boundary value prob- 
lems, in particular of the Dirichlet and Neumann problems, and give a moment 
problem which is responsible for this properties. Thus we will obtain some prob- 
lems that are equivalent to the problem (jl.ip . (|1.3p in the setting (|2.2|) . For details 
and generalizations see |16j or [17] . This equivalence is based on a connection 
condition of solution traces for the Cauchy problem. 

Let us consider a hyperbolic equation in an arbitrary bounded domain C M. 2 
with smooth boundary 

Lu = au'^ Xl + 6<x 2 + c< 2X2 = 0. (2.8) 

in the Sobolev space H m (il,),m > 3. 

Introduce a conormal vector and a derivative with respect to the conormal 
by means of an analog of the Green formula for the Laplace operator 

/ (Lu ■ v — u ■ Lv) dx — i (u' u Jv — uv'^) ds. 
Jn Jan 

One can count up that ^ = J(i/)£ - ^[l(v(s))]' s ■ §- a , 1(0 = a£ 2 + + c£f is 
the symbol of the operator L, v is an unit vector of normal, s is natural parameter 
on d£l, k = ±\v' s \ is the curvature, more exactly v' s = kr where r = (—1/2, v\) is the 
tangent vector. 

And let's consider an overdetermined boundary value problem for the equation 

u' s \aa = 7, u 'v, I an = «• (2.9) 
that is otherwise written down the Cauchy problem u\an = Jo, u'Jan = ko- 

We ask a question: what is a connection between the functions 7 and k if they 
are generated by a solution u of the equation (|2.8p ? In order to answer we need 
the following construction. 

The equation l|2.8p we write down also as 

(V ■ o 1 )^ ■ a 2 )u = (2.10) 

where a J = (a^a^^j = 1,2 are unit real vectors. Let's enter vectors a 1 = 
(—a\,a\), a 2 = (—a 2 ., a 2 ), directing vectors of set of characteristic directions A = 
A 1 U A 2 , AJ = {AfiJ|A € K}, j = 1,2, (a J ',a J ') = 0. Let ip Q = l Pi - ip 2 where (fj 
is any solution of equation tampj — Xj, i.e. tpj is an inclination angle of vector 
of the characteristic direction corresponding to the root Aj, ipo is angle between 
characteristics, and let A = sin^o = det jja 1 a 2 ||, aj are columns. Here and below 
a vector of a characteristic direction is understood as a vector v G C 2 which is a 
null of the symbol: l(y) = 0. The traces 7 and k of a solution u are linked by the 
following relations. 

Statement 1. If a function u E H m (Vt), m > 3 is a solution of the problem (|2.9[) 
for the equation IJ3HJ then the functions 7 € H m - 3 ^ 2 (dn), k G H m - 3 / 2 (dn) from 
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()2.9j) satisfy the conditions 

VQ GR[z], 



m 



<><> 



Q(x(s) • a 1 )^ = 0, 



Q{x(s) ■ a 2 )ds = 0, 



(2.11) 
(2.12) 



k e H m - 3 / 2 (dn),m > 3 satisfy 



VQ g M[ 
where x(s) is a moving point on 

An inverse statement holds also: 
Statement 2. If functions 7 G H m -^ 2 (dfl), 

the conditions (|2.11l) . (|2.12l) then there exists a solution it e iI m ~ 1_e (Q) of the 
problem (|2.9j) for the equation (|2.8p with each e > 0. For functions rp = u\gn, 
X = u'^ga we have l{v)x = K + [K^)]t/^^ an< ^ a ^ so ^ = / "f(s)ds + const (the 
Luzin's trigonometrical integral). In addition, function it is restored univalently up 
to an additive constant. The mapping: 

H m (dn) x ir"(<9ft)/{const} 9 {( 7) k) with (j2TT|l . (l2~T2ll } -> u G iT"- 1 ^^) 

is continuous (for all e > 0). 

Corollary 1. For each solution u G H m (f!,),m > 3 of the equations (|2.8[) the fol- 
lowing Zhukovsky's equality 



nds = 



(2.13) 



Oil 



holds. 



Proof. It follows from the condition (|2.11|1 for Q = 1 because j gQ r y(s)ds = 0. □ 



Consider following a moment problem: 

/ a(s)(x(s) ■ a j ) N ds = j = 1,2; N e Z+, 
Jan 

where on two given vectors a? G R 2 and on two sequences of numbers fj, 3 N it is 
found the function a. Obviously, for the case when dQ is the unit circle and 
vectors a? , j — 1,2 are equal a 1 = (1, i); a 2 = (1, — i) this moment problem turn on 
well-known trigonometric moment problem because then (x(s) ■ a J ) — exp(±iiV). 
Another way to the same is to write the Chebyshev polynomial Tjv insteed of the 
power. 

Among a lot of problems connected with above moment problem we will consider 
the problem of indeterminacy (uniqueness): for what curve dQ and vectors a J , j = 
1 , 2 there exists a function a such that 

V7V G Z+, j = 1,2, [ a(s)(x(s) ■ d r ) N ds = 0. (2.14) 
Jon 

The following fact takes place 

Statement 3. Let m > k > 3 and let we have three sets of statements: 

l m ) The homogeneous moment problem (|2.14j) has a nontrivial solution a G 
H m - 3 / 2 (dSl). 

2k) The Dirichlet problem u\dn = for the equation (|2.8p has a nontrivial 
solution u G H k (fl). 
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3/c) The Neumann problem vl v \on — for the equation (|2.8|) has a nonconstant 
solution u E H k (tt). 

Then l m ) 2 m _ g ); l m ) =*> 3 m _ g ); 2 m ) l m ); 3 m ) l m ) with <? = 1 + 
(By definition, for bounded domain H k+0 (tt) = \J f >0 H k+e (n)). 

Proof. 1) => 2). Using pair 7 = 0, k = 2a/ A with the help of the statement 1 we 
build the solution u € H m -i(dtt). 

2) => 1). We put a — k and apply the statement [2j 

The implications 1) => 3) and 3) 1) are similar. □ 

Note that instead of the Neumann problem in the statement 3 we can write the 
boundary condition of view 

u Vf = Aii 7 (2-15) 

with an arbitrary constant A. 

Note also that the case when the domain f2 is an ellipse has examined in the works 
by A. Huber [32], R. Alexandrjan [4j, V.I. Arnold 6! and one of the authors [16] 
or [17] . An answer to the question about properties of such the Dirichlet problem 
(|l.ip . (|1.3j) is the following. Reduce by means of a linear transform our problem in 
an ellipse to the problem (|1.3p in the unit disk for the equation (|2.10p . Find slope 
angles ip±, if 2 of characteristics and an angle tpo = — f2 between them. 

Statement 4. ([16], see also section [6} The problem (|l.lj) . (|1.3p has a nontrivial 
solution in the unit disk in a space H k (fl), k > 2 if and only if 

Vo/tt e Q. (2.16) 

If the condition (|2.16p is fulfilled then there is a denumerable set of linear indepen- 
dent polynomial solutions of the problem (|l.ip . (|1.3p . 

Thus, the existence of a nontrivial solution of the Dirichlet problem (|l.l|l . (|1.2|) in 
a general bounded domain with smooth boundary is equivalent to the existence of 
a nonconstant solution of the boundary value problem (|l.ip . (|2.15p . in particular of 
the Neumann problem, and is equivalent to the existence of a nontrivial solution of 
the moment problem (|2.14p . Below we will give a criterion of nontrivial solvability 
each of these problems with the curve (|1.4p in a view of the condition (|2.16p . 

2.4. Duality equation-domain and one more equivalent problem. Let ft C 

R™ be a bounded semialgebraic domain given by means of the inequality £1 = {x S 
W 1 : P(x) > 0} with a real polynomial P. Equation P(x) = gives us the boundary 
dtt. We assume the boundary of domain £1 is nondegenerate: | VP |^ on d£l. 
Consider the Dirichlet boundary value problem for the equation (|2.8p of the order 
2 with constant complex coefficients: 

Lu = L{D x )u(x) =0, u \ da = 0, (2.17) 

where D x — —id/dx. We understand a duality equation-domain as a correspon- 
dence the problem (|2.17p and the equation: 

P(-Dt){L(Ow(0} = 0, (2.18) 

given in the following statement: 

Statement 5. For each nontrivial solution of the problem (|2.17p from C 2 (f2) there 
exists unique nontrivial analytic in C solution of the equation (|2. 18|) from some 
class Z of entire functions and conversely: for each nonzero solution w € Z of the 
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equation (|2.18j) there exists a nonzero solution u G C 2 (£l) of the problem (|2.17p . 
The class Z is determined here as space of Fourier images of functions of the form 
9nV, where v G C 2 (K n ) and On is the characteristic function of the domain fl. 

For a clearness we give here a sketch of the proof. We assume the problem has 
a nontrivial solution u in C 2 (£l), let u G C 2 (R n ) be any smooth continuation of u 
on the R™ then we multiple u on characteristic function 8q of domain (On = 1 
in domain Q and Oq = out f2) and apply the operator L to the product Qq,u. 
Differentiating the product by Leibniz we obtain 

L(O a Z) = OnL(u) + L^u, Vu)S 8n + L 2 (u)(d d n)' v (2.19) 

where Sgn is measure supported on dfl: < 5qu,4> >= fgn^^^^' ^(u) = L(v)u, 
v is the external normal as earlier. The first term in (|2.19[) is equal to zero by 
means of equation. Take into account the boundary condition ugn = 0, obtain the 
last term will be transformed into term of view of the second term: L\(u 1 Vti) 6qq 
as, for example, for the case of one variable xS'(x) — —S. Then we multiple the 
obtained equality on P(x), the right-side part vanishes because e.g. for the case of 
one variable xS(x) — 0. We obtain an equation P(x)L(D)(0qu) — 0. Applying the 
Fourier transform T we obtain the equation (|2.18[) where w — J-(0q,v). Necessity 
is proved. Sufficiency will be obtained by means of conversion of this proof. A full 
proof for a general case see in the works [17] or [15] . 

The term 'duality equation- domain' means here an equivalence of the problems 
(|2~T7)1 and (j2~Tg]l that we will read here as: 

L(D x )u = 0, u \ P(X)=0 = 0, (2.20) 

P(-D 6 )v = 0, v |l(o=o= 0. (2.21) 
Now from statement [5] it follows 

Proposition 3. If there exists a nontrivial solution of the problem (|2.20|) with the 
string equation L(£) = — where £ = (C 1 ,? 2 ) is a covector, in a plane 

bounded domain with the biquadratic curve (jl.4p as a boundary dCl, then there 
exists also a nontrivial solution v G Z of the problem 

E a ^Sfe = °> u le=o = 0, v\ v=0 = (2.22) 

i,k=0 ^ ' 

and it is conversely. 

The last boundary value problem for the equation of the fourth order has only 
two boundary condition insteed of fourth condition as it should be e.g. for the 
problem of a Goursat (or Darboux) type. Therefore it is no wonder that a non- 
trivial solution of the problem (|2.22p exists there. But as we will show below for 
almost each curve (|1.4I) the problem Ijl.ip . (|1.2p has only trivial solution, therefore 
almost each problem (|2.22[) has only trivial solution as well. But this statement 
can seem astonishing namely because of the insufficiency of the data, in spite of the 
stipulation that the solution v belongs to the space Z. 

We finished the statement of propositions on boundary value problems for general 
domains and now we should wait for section 6 in order to obtain explicit answers 
for the Dirichlet problem in domains with biquadratic boundary. 
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3. Generic biquadratic curve 

3.1. Parameterizations of generic biquadratic curve. The complex curve 
(II. 4p is remarkable by that this is the most general algebraic curve with the prop- 
erty that almost each vertical or horizontal line intersecting C intersects it in 2 
points. 

Let (|1.4j) be a generic nondegenerate real biquadratic curve. Assume that the 
parameters are chosen such that the real curve C bounding the domain satisfies 
the condition (|2.ip of the subsection 12.21 

We will begin our study of the problem (|l.ip . (|1.3p by an observation that the 
curve (|1.4p is elliptic curve allowing an uniformization in terms of elliptic functions. 

Indeed, rewrite the equality (|1.4p in one of the form 

A 2 {x)y 2 + A 1 (x)y + A Q (x) = (3.1) 

or 

B 2 (y)x 2 + B 1 (y)x + B (y)=0, (3.2) 

where Ai(x), Bi(y) are polynomials of degrees < 2. We multiply the equation (|3.ip 
by A2, the equation (13. 2p by B2 and reduce these expressions to the forms 

Y 2 -D 1 (x)=0 or X 2 - D 2 (y) = 0, (3.3) 

where Y = 2A 2 {x)y+A 1 (x), X = 2B 2 (y)x+B 1 (y); {x,y) - (x,Y), (x,y) - (X,y) 
are birational transformations and Di(x), D 2 {y) are the discriminants of quadratic 
equations (|3.ip and (|3.2p : 

D x (x) = A 2 (x) - AA (x)A 2 (x), D 2 (y) = B 2 (y) - AB (y)B 2 (y). 

In a general situation the discriminants D\ 2 are polynomials of 4-th or 3-d degree. 
Recall [60] , that every curve of the kind 

Y 2 = iTiix) (3.4) 

with a generic 4-th degree polynomial n^x), can be transformed to a canonical 
form 

Y 2 = 4x 3 - g 2 x - 53 (3.5) 
which admits a standard parameterization 

x = p(t), Y = p'(t) (3.6) 

in terms of the elliptic Weierstrass function p(t) with primitive periods 2u)\, 2lo 2 . 
The parameters 32, 33 are so-called (relative) invariants of the polynomial D\. They 
are real for real 714. 

As {x, Y) — > {X, y) is a birational transformation, the primitive periods of the 
curves (|3.3p are coincide. The invariants 32,33 may be found by periods. Hence we 
obtain the following important statement (which was mentioned by Halphen [31]): 

Statement 6. Invariants g 2 ,g3 of polynomials D\(x) and D 2 (x) are the same. 

Thus, both curves (13. 3p are elliptic ones i 6Q] with appropriate primitive periods 
2u>x,2u2, the curve (|1.4p is homeomorphic to the torus: C w C/(2wiZ © 2^2^) and 
there are standard structures of a commutative group and an Abelian variety [29j . 
We deal with elliptic functions of the second order. Recall [60] that properties of 
a general elliptic function can be characterized by the number of poles (or that is 
equivalently, zeroes) in the fundamental parallelogram of periods. This number is 
called an order of the elliptic function (the order is taken with account of multiplicity 
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of poles). By the Liouville theorem, the simplest possible order is 2 |60J. E.g. the 
Weierstrass function pit) is of degree 2 because it has the only a pole of multiplicity 
2 at the point t = of the fundamental parallelogram. 

The general elliptic function of the second order has two arbitrary poles 
Pi , p 2 and two arbitrary zeroes Ci j C2 m the parallelogram of periods. The only 
condition is [60j p\ + p 2 — £i — ( 2 = 0, where f2 = 1m\L0\ + 2m 2 u 2 is an arbitrary 
period. It can be easily showed that generic elliptic function of the second order 
with given periods 2wi , 2^2 can be presented as 

= ap(t-t )+p 

Thus, $(t) depends on 4 independent parameters, say a, (3, 5, to. 

There is another, sometimes more convenient, representation of the function 

*(*): 

<r[t — di)a(t — d 2 ) 
where k is a constant and parameters e±, &i, d±, d 2 are related as 

e 1 +e 2 = d 1 +d 2 . (3.9) 

The form (|3 . 8[) is obtained from standard representation of the arbitrary elliptic 



function in terms of the Weierstrass sigma- functions [60] . In this case the points 
e\, e 2 and d±, d 2 coincide with zeros and poles of the function $(<) and the condition 
(|3.9p is equivalent to a balance condition between poles and zeros of the generic 
elliptic function. 

Note that apart from p(t) there are another special cases of functions of second 
order, e.g. the Jacobi elliptic functions sn(t; fc), cn(t; k), dn(t; k) [60j that we will 
use below. 

What is an uniformization for the biquadratic curve? The answer is given by 
following 

Theorem 1. The generic complex biquadratic curve (|1.4|) can be parameterized by 
a pair of elliptic functions of the second order with the same periods: 

x(t) = $!(*), y(t) = * 2 (t) (3.10) 

Conversely, any two elliptic functions x = 3>i(t), y = & 2 (t) of the second order with 



the same periods satisfy an equation ijl.4|) . 

This theorem is proved essentially e.g. in the famous Halphen monograph |31j 
on elliptic functions. We give a proof based on some main Halphen's ideas here. 

Proof. From (| 1 .4|) it follows the Euler differential equation ([33]) 

dx ± 



because dy/dx = -F x /F v = -X/Y (see ([33]) ). 

As we already saw, the polynomials D\(x) and D 2 {y) have the same invariants 
<72)<?3- Hence they both can be reduced to the same canonical Weierstrass form 
(see, e.g. [T], s.34) by means of a pair of the Mobius transforms 

x=f^, y = f^l, ^-^ = 1,1=1,2 (3.12) 
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with some complex parameters /ii, ... 772- Hence the equation (|3.11|) becomes 

dx dy 



yj 4£ 3 - g 2 x - g 3 \/ 4y 3 - g 2 y - g 3 



(3.13) 



because dy/dx = — X/Y as above. But the equation (|3.13|) means that for appro- 
priate periods 2us\ , 2u>i 

x = p(u), y = p(u + u ), (3.14) 

where p(u) = p(u; uj\, oj 2 ), u is an uniformization parameter and uq is a complex 
constant. Now we can return to initial variables x, y by means of inverse Mobius 
transforms to get finally 

^ = aip(u) + fa , y _ a 2 p(u + u ) + p 2 _ 15 > 

7ip(w) + <5i' 72p(u + u Q ) + 5 2 ' 

with some constants a.;, . . . , Si % . onbi — Piji = 1, i = 1,2. Then we apply formula 
(|3.7[) which says that we indeed obtained a pair of elliptic functions of the second 
order. 

The inverse statement of the theorem follows from a general theorem that any two 
elliptic functions x(t) and y(t) with the same periods satisfy an algebraic equation 
F(x(t), y{t)) = 0. The degrees of the polynomial F(x, y) with respect to variables x 
and y are determined by orders of corresponding elliptic functions. If both functions 
have the order 2 then the polynomial F[x, y) has the degree at most 2 with respect 
to each variable (this statement is contained in [60] as a problem for the reader). 
□ 

Thus, we proved the theorem [TJ Moreover, simultaneously we have proved the 
following 

Proposition 4. There exist complex linear- fractional transforms (|3.12[) such that 
transformed generic curve ljl.4p having the same form can be parameterized by 
only the Weierstrass function p(u) as in (|3. 14|) . 

Considering the formulae (|3.10[) we can write down our parameterization in terms 
of the Weierstrass sigma-function: 

X{t) = Kl a(t-di)a(t-d 2 y y{t) = K2 o-(t-d 3 Mt- di ) (3 ' 16) 

with two restrictions ei + e 2 = di + d 2 and 63 + 64 = ^3+ d±. If now we make a 
shift t — > t — to of the uniforming parameter t then we can choose to such that, say, 
e 2 = —ei,d 2 = —di. This means 

Proposition 5. In formulae (|3.10|l the function $i(t) can be chosen even: $i(— t) = 
<I>i (£) by means of choice a shift t — * t — to. 

From our considerations it follows an important corollary (which was mentioned 
by Halphen [3T] as well): 

Proposition 6. Consider the differential equation (|3.1ip . Let D\(x) and D 2 (y) be 
polynomials of degree 4 or 3 with the same invariants g 2 ,g 3 . And let (x(t),y(t)) 
be a solution of this equation (parameterized e.g. by an initial condition). Then 
x(t),y(t) satisfy a biquadratic equation of the form (|1.4[) . 
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3.2. Biquadratic foliation and singular points. There is an interesting me- 
chanical interpretation of last results. Assume that we have a dynamical Hamil- 
tonian system for two canonical variables x(t), y(t) satisfying a system of equations 

dH{x,y) . dH{x,y) 
X =—dy—> V = to - ' (3 ' 17) 

where H (x, y) is a Hamilton function of the system. Obviously H{x 1 y) is an integral 
of the system (13.171) . i.e. (H(x, y))' = 0. Choose the Hamiltonian as the biquadratic 
function (|1.4[) : H(x,y) = F(x,y). We have F(x,y) — c with some constant c de- 
pending on initial conditions for x, y. This constant can be incorporated to the coef- 
ficient aoo, so we can write down F(x, y) — 0, where F(x, y) = F(x, y) — c is again a 
biquadratic curve (note that for F(x,y) the coefficients A 2 {x) , A\(x) , B 2 {y) , B\{x) 
remain the same whereas the coefficients Ao(y) and Bo(y) differ from initial ones 
by a constant). Then 

9H i X,V) = 2A 2 (x)y + A x {x) = ±JK(x), 
dy v 

where D±(x) — Af(x) — 4A 2 (x)Ao(x) (because y can be excluded as a root of 
quadratic equation A 2 (x)y 2 + Ai(x)y + A (x) = 0). Quite analogously 

dH(x,y) 



ih -2B 2 (y)x + B 1 (y)=±^D 2 (y). 

We thus see that for any fixed Hamiltonian level H = c the variables x(i), y(t) satisfy 
the Euler equation (|3 . 1 1 [) where polynomials Di(x), D 2 (y) have the same invariants 
!J2, ff3- Note, that in this case the invariants g 2 , <?3 (and hence the periods 2u)\, 2uj 2 ) 
will depend on the value of integral c. We thus obtain a whole one-parameter family 
of biquadratic curves F{x,y) = c and corresponding elliptic functions x(t),y(t) of 
the second order that are trajectories of this dynamical system. 
Rewrite discriminants Di(x), D 2 (y) in factorized forms 

4 4 

Dx{x) = qi Y[(x - Xi), D 2 (y) = q 2 Y[{y - Vi), 

i=l i=l 

where q±, q 2 are leading coefficients of the discriminants and x,, y^ i = 1, 2, 3,4 are 
their roots (for simplicity we assume that both discriminants have degree 4). In 
general, roots Xi,yi will depend on the parameter c. What is mechanical meaning 
of points Xi, From equations of motion it is seen that 



± 



\Jdi(x), y = ±^D 2 (y) 



Thus Xi and yi are stable points: x\ — jji — 0. We should demand that the points 
(xi,yk) belong to our biquadratic curve F(x,y) = 0. This leads to conditions that 
points (xi,yk) satisfy the conditions 

F(x, y) = 0, d x F{x, y) = 0, d y F{x, y) = (3.18) 

which in turn means that points (xi, yk) are singular points of the biquadratic curve 
F{x,y) = 0. As we saw, in generic situation this curve is elliptic and hence has a 
genus 1 (genus C = C^_ 1 — d, n is degree and d is a number of double points of 
C). The degree of this curve is 4. Assume that the curve is irreducible. Then such 
curve cannot have more than 3 singular points in complex projective plane. The 
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latter is defined by the coordinates (so,si,S2) such that x — s\/sQ,y — s 2 /sq. In 
these coordinates we have equation of our curve 

a 22 s\sl + s sis 2 {ai 2 si + a 21 s 2 ) + sl(a 20 s\ + a 02 s 2 2 + a n sis 2 ) + 
So(aioSi + aoiS2) + Sgaoo = 0. (3.f9) 

Elementary considerations show that two points (0, 1, 0) and (0, 0, 1) of the pro- 
jective plane are singular for any values of parameters a^. Thus, only one singular 
point can exist in each finite part of the plane. In turn, this means that at least 
two roots, say xi, x 2 of the discriminant D\(x) should coincide: x\ = x 2 . The same 
condition holds for the discriminant D 2 (y), i.e. yi = y 2 because invariants (72,33 
of the discriminants D\(x) and D 2 (y) are the same. Then it is elementary verified 
that the point (xi,y±) will be indeed a singular point of the biquadratic curve. In 
principle, the second singular point can occur. But in this case the genus of the 
curve will be -1 meaning that the curve is reducible. 

We thus have the following 

Proposition 7. The irreducible biquadratic curve F(x,y) = has the a singular 
point in a finite part of the (complex) plane if and only if the discriminant D\(x) 
(and hence D 2 (y) as well) has a multiple zero x\ (correspondingly yi). In this case 
the point (x\^y\) is singular and unique. 

Note this proposition can be reformulated in an equivalent form. Indeed, the 
polynomial D\(x) has a multiple zero if and only if its discriminant is zero. Thus, 
in order to find all singular points of the curve we should first calculate the discrimi- 
nant D\{x) of the biquadratic curve F (x, y) = and then calculate the discriminant 
A(Di(x)) from the polynomial Di(x) ("discriminant from the discriminant"). We 
have obviously coincidence of the two discriminants 

A(D 1 (x)) = A(D 2 (y)) 

because the invariants (72,53 of the polynomials Di(x) and D 2 (y) are the same 
(statement [6]). Condition A(Di(x)) = (or, equivalently, A(D 2 (y)) = 0) gives 
us some nonlinear equations for the coefficients ciik- Under such condition the 
biquadratic curve F{x,y) = has a genus < 1, i.e. it is either irreducible and has 
the only singular point in a finite part of the projective plane, or it is reducible: 
F(x,y) = Ti(x,y)r 2 (x,y), where Ti(x,y),T 2 (x,y) are two polynomials linear with 
respect to each argument x, y (but t^(x, y) arc not in general linear functions, they 
may contain the terms like xy). We have obtained 

Proposition 8. The condition A(Di(x)) (= A(D 2 (y))) ^ is necessary and suffi- 
cient for the equality: genus C = 1. 

3.3. Case of generic symmetric curve. Above we considered generic case when 
our biquadratic curve F(x 7 y) = is non-symmetric. 

Assume now that our curve is symmetric, i.e. F{x,y) — F{y 1 x). Equivalently, 
this means the coefficient matrix a%k in (jl.4p is symmetric: et^ = a^. In this case, 
obviously, both discriminants coincide D\ (x) — D 2 (x) . From the Euler differential 
equation (13. lip we obtain that a parameterization can be given by formulae 



x{u) = $(tt), y(u) = $(u + u ), 



(3.20) 
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where uq is a constant and $(u) is an even function of the second order, i.e. <I>(i) 
can be presented in the form 

jp(u) + S 

Thus, in the symmetric case a parameterization is provided by some even elliptic 
function of the second order. Vice versa, any pair (x, y) = ($(w), $(u + uq)) with 
arbitrary uq generates a symmetric biquadratic curve by formulas (|3.20p because 
the point (y, x) = + u a ), $(u)) = ($(■&), $({t + u )), u = —u — uq belongs C 
also. Note that the last statement was attributed to Euler in the work [55] , 
We obtain the following 

Proposition 9. The generic complex symmetric biquadratic curve (|1.4[) can be pa- 
rameterized by an even elliptic function of the second order and a shift as in (|3.20|) . 
Conversely, for each elliptic function (even or not) of the second order and any 
shift Mo the variables x, y from (|3.20p satisfy the equation (| 1 .4|) with a symmetric 
matrix A. 

This implies the important 

Proposition 10. The generic complex biquadratic curve (|1.4p can be transformed to 
symmetric one by means of a linear-fractional complex changes of variables (|3.12p . 
If the initial curve is real then it can be transformed to real symmetrical curve also 
(although corresponding transformation can be with complex coefficients). 

Proof. As we have shown in the thcorcm[T](see proposition^), generic non-symmet- 
ric curve can be transformed to a curve described as (|3.14[) by linear-fractional 
complex changes (|3.14l) that means the curve-image is symmetric one in virtue 
of the last proposition. In the real case we note that as it is well-known, the 
invariants 32,33 of real polynomial D\(x) are real, so that the differential equation 
()3.13|) is real and there exists its real solution (x(u),y(u)) that can be extended 
onto complex domain and therefore satisfies the symmetric equation (|1.4p . The 
last equation has irreducible polynomial which must be real if 022 = 1 because 
we always can choise real parameters u±, ...ug such that vectors Vi composed of 
components x k (ui)y l (ui), < k,l < 2,0 < k + I < 4 will be linear independent, 
then the coefficients of our polynomial satisfy a linear system of 8 linear equations 
with real coefficients (Vi)j and the real right-side parts —x 2 (ui)y 2 (ui). □ 

Note that another proof of the last fact is there in the work [35] , see below statement 

ra 

Let us now consider the question how restore the polynomial F if we know the 
discriminants D\,D2- In symmetric case T. Stieltjes 55 proposed a nice explicit 
formula for the polynomials F(x, y) by means of a solution of the differential equa- 
tion (|3.1ip . Assume that 

D 1 (x) = D 2 {x) = + 4&ix 3 + 6b 2 x 2 + 4b 3 x + b 4 . 

Then F(x,y) can be given by the determinant: 

1 -{x + y)/2 xy 

1 b 61 6 2 - 2C 
-{x + y)/2 hi b 2 +C 63 

xy 62 — 2C 63 64 



F(x,y) 



(3.22) 
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The function F(x,y) is denned up to an arbitrary non zero number factor. It is 
assumed that the curve is nondegenerated, i.e. its genus is 1. This is possible if 
and only if the determinant 

6 61 b 2 - 2C 

61 b 2 + C 63 
b 2 - 2C 63 b 4 

is nonzero. C is an integration (arbitrary) constant (|55j). The Stieltjes formula is 
useful in the problem of reducing of the arbitrary symmetric biquadratic curve to 
some standard forms, as we will see below. 

3.4. Canonical forms of biquadratic curve. The general curve (|1.4| contains 
8 free parameters. It is naturally to transform this curve to the form containing 
the smallest possible numbers of parameters. 

First, note that under arbitrary projective transformations 

x^^ X + m , x^ ^ X + 112 (3.23) 
H\x + V\ \x 2 x + Vi 

with complex parameters we obtain similar equation (jl.4p but with transformed 
parameters a^. This idea has already exploited in the proof the theorem 2. Ev- 
ery projective transformation (|3.23[) contains 3 independent parameters, hence it is 
possible to reduce the total number of independent parameters to 8 — 6 = 2. As 
our curve is an elliptic one these free parameters are only invariants g 2 , 33 under 

linear-fractional transformations of variables x and y separately. More explicitly, 

2 

if we consider the projectivisation a i\i23\h x \ x 2 V 3 \ vl? ^ x — x iiV = Vi 

*ll*2 >jl ,32 = 1 

of initial curve l|1.4[) and its projective transformations of variables x and y sep- 
arately that <?2j<?3 are only invariants (of respective wights 4 and 6) but for the 
transformations of the group SX(2,C) these are absolute invariants. 

Above we reduce the general case to the symmetric one so that we restrict 
ourselves with considering a symmetric curve F(x,y) — F(y,x). So we would like 
to find some canonical forms of the curve containing only 2 parameters. 

There are two obvious canonical forms which can be obtained for the curve 
F(x, y) = 0. These two forms correspond to two canonical forms of elliptic integrals 
in the Euler differential equation (|3.11j) . 

(I) The first one can be obtained if one reduces polynomial D\(x) — D 2 {x) to 
the canonical cubic Weierstrass form: 

Di(x) = Ax 3 ~ g 2 x - g 3 . 

Such form can be always achieved by an appropriate projective transformation (with 
possible complex coefficients). Then from the Stieltjes formula (|3.22[) we obtain the 
expression (see also [31] ) 

F(x, y) = (xy + (x + y)w + g 2 /A) 2 - (x + y + w){Axyw - g 3 ) = 0, (3.24) 

where w = C is an arbitrary parameter. There is a parameterization of this curve 
in terms of the Weierstrass elliptic function 



x{u) = p(u), y{u) = p(u + u ) 
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and w = p(uo), where uo is an arbitrary constant. It is easy to calculate the 
discriminant: 

Dx{x) = (4iu 3 - g 2 w - g 3 ) (4x 3 - g 2 x - g 3 ) 

If the parameter w is such that 4ui 3 — g 2 w — 53 = then D\(x) = and in this case 
the curve F(x,y) is reducible: F(x,y) — pr(x,y), where p(x,y) is a polynomial of 
degree 1 with respect to each variable x, y. If 4w 3 — g 2 w — 53 7^ then the curve 
is irreducible. The singular points in a finite part of the complex projective plane 
appear only if g 2 = 11g\. This condition means that the polynomial 4x 3 — g 2 x—g 3 — 
A(x — ei){x — e 2 )(x — e^) has a double root, say ei = e 2 . In this case we can put 
g 2 = 3r 2 , ga = r 3 , where r is some parameter. Then it is elementary verified 
that the point x — y = —t/2 will be the only finite singular point of the curve 
F(x, y) = 0. If g 2 — 27gf 7^ then there are no singular points in the finite part of 
the projective plane and the curve (|3.24[) is irreducible and has a genus 1. 

The curve (|3.24|) contains 3 parameters w, g 2 ,gs- Assume that 4w 3 — g 2 w — 53 = 
W 7^ 0. In this case the curve is irreducible. By the linear transformation of 
arguments x — > ax + w, y — > ay + w, where a 3 = 1/W we can eliminate terms 
x 2 y, xy 2 and x 2 + y 2 . The curve then is reduced to the form 

F(x, y) = x 2 y 2 -x-y + Axy + B = (3.25) 

containing only 2 independent parameters A, B. The curve is elliptic non-singular 
(i.e. having genus 1) if condition 

A(Z>x) = B(A 2 - AB) 2 + A{36B - A 2 ) - 27 ^ (3.26) 

holds. 

Thus we proved the following 

Proposition 11. The generic complex biquadratic curve (jl.4|) can be transformed to 
canonical form (|3.25p by means of a linear-fractional complex changes of variables 
(|3.12jl . If the initial curve is real then it can be transformed to real form also 
(although corresponding transformation can be with complex coefficients). 

As we saw, any irreducible biquadratic curve can be transformed to this form. 
However, here we need general projective transformations (|3.23|) with complex pa- 
rameters. This will be so even in the case when all parameters of the biquadratic 
curve are real. 

(II) In beginning let remember well-known Legendre transformation. Assume 
that an elliptic curve has already the form Y 2 = ^4(2;) (see (|3.3li l. where ir^x) 
is an arbitrary polynomial of the 4-th degree with real coefficients. Using only 
a linear- fractional transformation x — T(x) with real coefficients it is possible to 
reduce this curve to the form (with a new Y) where tt4,(x) has only even degrees of 
x: 7r 4 (x) = ax 4 + fix 2 + 7 with some real coefficients a, f3, 7 (see, e.g. [2"3"]V 

Assume now that our biquadratic curve is symmetric and has only real coeffi- 
cients ciik- Then we can always transform this curve to a form without terms of 
odd degrees, i.e. 021 = a%2 = aio = aoi = 0. Indeed, taking x = y in expression 
(jl.4p for F we obtain the polynomial F(x, x) of fourth degree and apply to it the 
Legendre change x = T(x) in a generic case that reduces the equation F(x, x) = 
to the view ax 4 + (3x 2 +7 = 0. Therefore the change with the same V 



x = T(x), y = T(y) 



(3.27) 
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applied to (|1.4|) gives us an equation with desired polynomial (observation from 
[35]). 

We then obtain so called the Euler-Baxter biquadratic [TO], [55] : 

F(x, y) = x 2 y 2 + a(x 2 + y 2 ) + 2bxy + c = 0, (3.28) 

where a, 6, c are remaining real parameters if the initial curve was real and, obvious, 
will be complex if the initial curve was complex. The curve (|3.28|) plays the crucial 
role in deriving addition theorem for the elliptic function sn(t) pp. It appeared also 
in Baxter's approach to 8- vertex model in statistical mechanics [TP] . 

We first analyze possible finite singular points of the curve (|3.28|) . The discrim- 
inant of the curve (|3 . 28[) is equal to 

U2 _ 2 _ 

D^x) = ~ax 4 - (a 2 -b 2 + c)x 2 - ac = -a(x 4 - bx 2 + c); b = , (3.29) 

a 

Di(x) = D 2 (x) and by proposition [7] the curve (|3.28[) has a finite singular point 
under the condition 

A(Di(x)) = I6a 2 c {{a 2 - b 2 + cf - 4a 2 c) 2 = 0. 

If c = then we have an obvious singular point x = y = 0. This singular point 
will be isolated for all values of the parameters a, b excepting the case a ± b = 0. If 
a ± b = then the curve F = becomes reducible. In any case if c = the genus 
of the curve (|3.28p is less then 1 (i.e. for generic case and —1 for exceptional case 
corresponding to a reducible curve). Thus in this case the curve is not elliptic and 
can parameterized by rational functions. Here the change x — > k/x, y — > n/y gives 
us a form F = x 2 + y 2 + axy ± 1. 

Using then scaling transformations x — > kx, y — > ny in real case we can reduce 
the coefficients c to ±1 depending on sign of this coefficient. 

Thus, in real symmetric case we have simplest forms: 

(i) F = x 2 y 2 + a(x 2 + y 2 ) + 2bxy + 1; 

(ii) F = x 2 y 2 + a(x 2 + y 2 ) + 2bxy - 1. 

In (ii) one can consider a > by virtue of substitution x — > y — > 1/y. 

The canonical biquadratic curve (|1.4p with F from (i) or (ii) we will call the 
Euler-Baxter curve (as (|3.28p V 

Proposition 12. The generic complex biquadratic curve (|1.4[) can be transformed 
to canonical form (i) by means of a linear-fractional complex changes of variables 
(|3.12[) . If the initial curve is real then it can be transformed to real form (i) or 
(ii) (although corresponding transformation can be with complex coefficients). By 
means of only some real linear-fractional transformation (|3.12p any generic real 
symmetric biquadratic curve (|1.4|) can be reduced to one of two form: (i), (ii). 

Remark that in the paper |35] it is proved the following 

Statement 7. By means of only some real birational transformation any generic real 
elliptic biquadratic curve (jl.4p can be reduced to one of three form: (i), (ii) or to 
the form 

(iii) F = x 2 y 2 + a(x 2 - y 2 ) + 2bxy - 1. 

Note that the discriminants of the curve (iii) are equal to 

DAx) = -a(x 4 - l>x 2 + 1); D 2 (y) = a(y 4 + by 2 + 1): b = - + ^ + 1 . (3.30) 

a 
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In the complex case we can reduce the coefficients c 7^ to 1 and leave only the 
case (i) among nondegenerate cases. Here following |10] we can find a parameteri- 
zation of the curve (|3.28p in terms of elliptic Jacobi function sn: 

Proposition 13. The curve Q3.28P can be parameterized by the formulae 

x = V~ksn(t;k), y = Vk sn(t ± 77; k) , (3.31) 

with any sigh ±, where parameters k,r] are determined by relations 

k + k- 1 = (b 2 -a 2 -l)/a, 1 + kasn 2 (r ri k) = 0. (3.32) 

Note that we can replace t = t + K in order to deal with an even function. 
Note that the mapping (|3.3ip is an analytic diffeomorphism from fundamental 
parallelogram, factored by a standard way into torus, onto C. 

Real paramctcrizations of curves (i), (ii) and (iii) are brought in [35] and [36) . 

Now let's remember expressions p. 31) writing as 

F 2 ( Xl , yi ) = D 2 ( Vl ) = B 2 1 (y 1 )-AB (y 1 )B 2 (y 1 ). (3.33) 

We deduce y-extreme point (xi,j/i) gives us = y'(xi) — F x /F y (x\, yi); = 
Fx(xi,Vi) = 2B 2 (yi)x 1 +B 1 (y 1 ). If B 2 (yi) = then an = -B (l/i)/Bi(j/i) from 
(O, otherwise F 2 (x uVl ) = D 2 {yi) = 0, xi = -B 1 {y 1 )/2B 2 {y 1 ). Thus, for any 
y- vertex (xi,yi) the number y\ is a root of the polynomial D 2 {y\). Conversely, 
from the equalities (|3.33|) and y'F y — —F x we see that for each such root yi of 
D 2 the point {x\, y\) is cither y- vertex (extreme point in y-direction) or a singular 
point. But as we have seen in subsection 13.21 a singular point can be only if a root 
of D 2 is multiple, i.e. A = g 2 — 27 g 2 = 0; this is an exceptional case. Therefore it 
is proved the following 

Proposition 14. For generic case the roots of D 2 (y\) and only they are y-coordinates 
of y-vertexes of the curve. It's analogous, for generic case the roots of Di(x\) and 
only they are ^-coordinates of x-vertexes. 

Now we can give geometrical an interpretation of cases (i), (ii), (iii), a proof of 
which can be easy obtained by analysis of formulae (|3.29p and (|3.30| . 

Proposition 15. 0). The real curve in cases (i), (ii) and (iii) disappears only in the 
case (i), b < 2, a > that is equivalent to b 2 < (a + l) 2 , a > 0. 

1) . In case (i) of nonvanishing real curve the equation D\(x) = (and the 
equation D 2 (y) = 0) has no real root only in the case b < 2, a < (there is no 
vertex) . 

2) . In case (i) the equation Di(x) = (and the equation D 2 (y) = 0) has four 
real roots under b > 2 (there are 4 £- vertices and 4 y- vertices) . 

3) . In the case (ii) the equation D\(x) = (and the equation D 2 (y) = 0) has 
two real roots (there are 2 a;- vertices and 2 y- vertices) . 

4) . In the case (iii) for a > (b > 2) the equation Di(x) = has 4 real roots and 
the equation D 2 (y) = has no real root (there are 4 x- vertices and no y- vertex). 

5) . In the case (iii) for a < (b < —2) the equation D\(x) = has no real 
root and the equation D 2 (y) — has four real roots (there are 4 y- vertices and no 
x- vertex) . 

If we will allow values x — 00, y — 00 i.e. consider our real curve on a thorus 
S 1 x S 1 then the number of x-vertexes and the number of y-vertexes will not be 
changed by Mobius transformations (|3.23[) . therefore we obtain 
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Proposition 16. Linear-fractional changes (13.23|) don't remove the curve from its 
class 0)-5) of the proposition fT5l 

3.5. John mapping of biquadratic curve and periodicity. At first, we show 
how complex John algorithm works for the curve C. Consider a case of some more 
general curve. 

Let C C C 2 be a complex curve described parametrically as 

x(t) = 4>(t), y(t) = <t>{t + e), teC (3.34) 

where e is a nonzero complex parameter. Assume that 4>{z) is an even periodic 
meromorphic function, i.e. 4>{—z) = 4>{z), (f>(z + T) = 4>(z) with some (complex) 
constant T. And let C = C D K 2 be a real curve given by means of a contraction 
of the functions x, y on some set S £ C. 

We will assume that the curve C satisfies the condition: 

The curve C is nondegenerate and each straight line of the form x = xq 

or y — ?/o intersects the curve C at no more than two points (3.35) 

Replace t — » —t. Then, due to the fact that <j)(t) is even we see that x(— t) = x(t) 
but in general y(—t) ^ y(t). This means that the point (x(t),y(—t)) on the curve 
is obtained as the second intersection point of C with the "vertical" line passed 
through the initial point (x,y). Thus the transformation i — > — t is equivalent to 
the involution I\ in the John algorithm. Quite analogously, the transformation 
t — > — t — 2e remains the coordinate y on the curve C unchanged whereas x is 
transformed to another point on the intersection of the curve C with "horizontal" 
line. Thus, transformation t —> — t — 2e is equivalent to second involution I2 in the 
John algorithm. For the John mapping T we have obviously 

T = hh <-> * -> t - 2e, T- 1 =hh ^t^t + 2e, (3.36) 
Thus, we proved the following 

Proposition 17. For curve C given by (|3.34[) with an even periodic meromorphic 
function <\> satisfying the condition (|3.35|) the complex John mapping T is equivalent 
to a shift of the parameter t by the step — 2e. 

Now the periodicity condition of the complex John mapping is obtained the form: 

2ne = mT (3.37) 

with some integer n, m. 

In the case of the biquadratic curve we have two periods 2ui, 2wi and then the 
periodicity condition is taken the form: 

ne = m\U\ + TO2W2 (3.38) 

with some integer n, mi, m.2- 

Now we have only two cases: either all points of C are not periodic (for the John 
mapping) or each point is periodic with the same period (that is equal to n if the 
condition © or (|3"38) is fulfilled). 

Consider the last canonical form of previous subsection - the Euler-Baxter curve 
(|3.28p that is parameterized as (|3.31[) . It is important that for given x — \fk sn(t; k) 
there are two values of y corresponding to two values of rj: y\ = \f~k, sn{t + tj; k) and 
j/2 = \/ksn{t — rj; k). These two values 2/1,2 correspond to two points of intersection 
of the line x = xq with the curve (|3.28p . 
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Thus we can find all points M n of the complex John algorithm: 

M n = (Vk sn(t + 2-qn - 3r/; k), \fk sn(t + 2<qn - 2-q; k)) (3.39) 
The periodicity condition with a period n is 

2ry n = AKmi + 2iK'm 2 (3.40) 



with some integer rai,m2 and constants K = K(k) = J dtj — k 2 t 2 )(l — t 2 ), 

o 

K' = K(k'), k' 2 = 1 — k 2 by well-known notations. 

In order to write analogous formulae in real case we should write real paramer- 
izations by means of real-valued elliptic functions. 

Consider first the case a > 0, c = +1 of real x-y— symmetric curve (|3.28[) 
and use the parameterization (|3.31|) . Above, in the proposition [TBI (s. 0) we noted 
that under a > the condition b 2 > (a + 1) 2 (or b > 2) is necessary and sufficient for 
the existence of real curve (|3.31[) . Therefore assume b 2 > (a + l) 2 then k + k^ 1 > 2 
and we may choose k = k\, < k\ < 1. From the second equality (I3.32|) we 
obtain that sn(r], k) is pure imaginary. This means rj = 2mK + 9i with an integer 
m and real 6 (see [5]). Moreover, x must be real, that is either t = nK'i + r or 
t = (2n + \)K + it with an integer n and real r ([9]). But y must also be real, 
i.e. t + 77 = K'n\i +t\ or t + r/ = (2n\ + 1)K + It\ with an integer n\ and real t\. 
Let us add expressions for t and 77, obtain that it is possible only a variant with 
the second formulae for t and t + r\: t\ = 9 + r. Thus we have suppose that the 
parameter t isn't real: t = ±K + zt, however the values x, y in the paramentization 
(13.31)1 are real and the parameter r is real, here the sign before K determines a 
branch of the curve which is two bounded ovals in this case. The relations (I3.32p 
have a denumerable set of solutions "q, we choose 9 = Im 77 as a minimal positive 
number of all Im 77. Now the periodicity condition (|3.40() can be written in the form 
(to = TO2): 

4 = -eQ (3.41) 

K' n 

where from relations (|3.32[) we obtain that the number 9 is a minimal positive 
solution of the equation sc(9,k') = \j\fak (here and below sc — sn/cn, ns — 
l/sn and so on as usually). We proved the foolowing 

Proposition 18. The condition (|3.41[) is a criterion of periodicity of the John map- 
ping T in the case a > 0, c = +1 of real x-y— symmetric curve p.28[) . Under it all 
points of our curve are periodic with the same minimal period n in nonreducible 
fraction — . The number to is a number of full turns of the curve that the mapping 
T n do. 

Consider also another cases of real x-y— symmetric curves (|3.28|1 . We take 
advantage of a list of cases given in the work [35] (note that the last work constains 
mistakes in above case c = 1, a > 0, here we used formulae from [10)). We write 
final parameterization formulae where the modulus k is given by b (k 12 + k 2 = 1)) 
and the shift 77 — by a or b (they are coordinated). Here b = (b 2 — a 2 — c)/a. 



Case a > 0, c = — 1, : x = yk/k' cn(t;k), y = \/kjk! cn(t ± 77; k), (342) 
where k/k' — k'/k = b, a ~ ds 2 (rj,k)/kk' , b=—cs(rj,k)ns(r],k)/kk'. 
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Case a < 0, c = — 1, : x = y/k/k' nc(t; k), y — y/k/k' nc(t ± 77; k), 
where k'/k—k/k' = b, a = —ds{r],k)/kk', b = cs(rj,k)ns(r],k)/kk' . 



Case a < 0, c = 1, b < —2 : x = y/l/k'cs(t; k), y = y/l/k' cs(t ± 77; k), 
where —l/k' — k' = b, a — ~cs 2 (i],k)/k' , b = ds(rj, k) ns(rj, k)/k' . 

Case a<0,c=l, b>2: 

unbounded part x = y/l/k ns(t; k), y = y/l/k ns(t ± r/; k), 
bounded part x — Vk sn(t; k), y = Vk sn{t ± 77; k), 

where 1/k + k = 0, a = —ns 2 (r],k)/k, b = cs(r),k) ds(r),k)/k. 



(3.43) 
(3.44) 



(3.45) 



Separation in the last case is connected with that the curve has a bounded oval 
and other unbounded branches. Last four cases are given in the work [35] - 

Remaining cases are considered in the work [36] : Case a < 0, c = 1, |6| < 2 of 
real symmetric curves (|3.28|) can be transformed to the case (|3.44p by means of 
substitution 

fl-x l-y\ m1 l + a-b - 2(l-o) ? „ 16a 
(x,y)= — — ,— . Then,a= — -=-, 6=-^ r4= & = 2J 



1+x'l + y/' ' l + a + &' l + a + 6' & 2 - (a + l) 2 



Case a < 0, c = 1, \b\ < 2 : x — Vk/ sett; k), y = y/l/k' sc(t ± 77; fc), 

(3-46) 

where — l/k' — k' = b, a = —cs 2 (r],k)/k' , b = ds(i~i,k)ns(r],k)/k' . 
For real x-y— asymmetric curves (|1.4p . the case (iii) (see the statement [7]) 



Case o > 0, c = 1, (6 > 2) : z = Vk'nd{t; k), y = ± vT/fc 7 dn(i ± 77; fc), ^ ^ 
where l/fc + fc = o, a = k'nd 2 {rj,k) 1 b = k 2 sd(rj,k) cd(r],k). 



Case a < 0, c = 1, (6 < -2) : x = Vk' sc(t; k), y = y/l/k' cs(±t + 77; k), ^ 

where — 1/k — k = b, a = —k' nd 2 (r], k), b = k sd(rj,k) cd(rj,k). 

All these cases (|3.42[) - (|3.48[) admit calculations as above that gave us the formula 
(|3.4ip . Such calculations were done by our collaborators M.V. Belogljadov and A. A. 
Telitsyna ([H],[5S]) and they bring the following results: 

a common answer for cases (|3~4^t . ([3~l4]) . ([3~4"5]) . (|3~^5|) . (|3~4^|) 



Re 77 m _ ,„ 

"2/T - n ^ Q; (3 - 49) 

= -eQ; (3.50) 

Im 77 , m „ , 

^/-l2 = -eQ. 3.51 
2iv ' 71 

Here fc and 77 are counted up from relations (|3.42[) - (|3.48[) corresponding to con- 
sidered case. Thus, the following proposition takes place: 



for the case (|3.43p 
for the case p. 471) 
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Proposition 19. The conditions (j3.49[) - (j3.51|) are criterions of periodicity of the John 
mapping T in corresponding real cases (|3.42|) - (|3.48|) . The number n is the period 
of dynamical system and m is the number of full turns of mapping T n . 

Note that a correspondence to cases of proposition [15] can be easy seen and 
presence of several cases corresponding to a case of the proposition [TBI is explained 
by technical reasons (by choice of an infinite point on projective line). 



4. The Poncelet problem 

In this section we demonstrate a nice correspondence between the mapping by 
F.John and famous the Poncelet problem (the Poncelet porism) for two conies. We 
start by recalling the Poncelet porism in a well-known form. 

4.1. The Poncelet porism in form of two circles. At the beginning let a circle 
A lie inside another circle B. From any point on B, draw a tangent to A and extend 
it to B. From the intersection point, draw another tangent, etc. For n tangents, 
the result is called an n-sided Poncelet transverse. This Poncelet transverse can be 
closed for one point of origin, i.e. there exists one circuminscribed (simultaneously 
inscribed in the outer and circumscribed on the inner) n-gon. We could begin with 
a polygon that is understood as the union of a set of straight lines sequentially 
jointing a given cyclic sequence of points (vertices) on the plane. If there exist two 
circles, inscribed and circumscribed for this polygon, then this polygon is called a 
bicentric polygon. Note that sides of the polygon can intersect and the intersection 
point is not obligatory to be a vertex. Furthermore, the inscribed circle does not 
obligatory touch a segment between vertices, the contact point can lie on extension 
of the side and therefore the circles can intersect. Bicentric polygons are popular 
objects of investigations in geometry. This is most known form of the Poncelet 
porism. If we denote by r the radius of the inscribed circle, by R the radius of the 
circumscribed circle and by d a distance between the circumcenter and incenter for a 
bicentric polygon then these three numbers can not be arbitrary and together with 
n they satisfy some relations. So, for the case of triangle the relation is sometimes 
known as the Euler triangle formula R 2 — 2Rr — d 2 = 0. One of popular notations 
for such relations (necessary and sufficient for existence of a bicentric polygon) is 
given in terms of additional quantities 

1 u 1 1 
b = — : , c = - . 



R+d ' R-d r 

So, for a triangle above the Euler formula has the view: a + b = c, for a bicentric 
quadrilateral, the radii and distance are connected by the equation a 2 + b 2 = c? . 
The relationship for a bicentric pentagon is 4(a 3 + b 3 + c 3 ) = (a + b + c) 3 . In a 
general case one introduses numbers 

2n 2 (a 2 - b 2 ) 

A = 1 + —rhi or , w = costr 1 A , k 2 = 1 - e - 2u , K = KVk) (see s. 

a z (b z — c l ) 

and then the relationship can be written by means of elliptic functions in the form 



K ,\ c^b 2 -a 2 + bV' 



.2 



sc - fc =- irjri ( 4>1 ) 

(Richelot (1830) - the first edition of the criterion, Kerawala (1947) ([37]) - the 
written criterion). 
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4.2. Setting of the Poncelet problem. Recall the Poncelet problem [13] for 
the case of two ellipses, for simplicity and as it was introduced by Jean- Victor 
Poncelet himself. We take two arbitrary ellipses A and B, A inside B. Let us 
have an arbitrary point Q\ on the ellipse A and pass a tangent straight line to A 
at the point Q\, This tangent crosses the ellipse B at two points Pi and Pi, Pi 
before Pi with respect to a standard orientation. Then we take the point Pi on 
B and pass the second tangent to the ellipse A. We denote as Q2 the point on 
A where this tangent contacts with A. This tangent meets the ellipse B in two 
points Pi and P3. Take the point P3 and repeat this procedure. Then we obtain a 
mapping Ub ■ B — ► B which acts by the rule Ub ■ Pfe — ► Pfe+i that will be called 
the Poncelet mappings below. Moreover, we obtain the mapping Ua ■ A — ► A 
acting by the rule Ua ■ Qk — * Qk+i- More precisely, because the definition of 
Ub exploits a point Q € A we should introduce two mappings I a, 1b ■ C — * C, 
C := {(Q, P) £ i x B I P lies on tangent line to A at Q}, acting by the rules 1a '■ 
(Qij-Pa) -> (Qa,^), : (Qi,-Pi) -* (Qi,p2)- These mappings 1 b Ja generate a 
composition U = Is Ia that is similar to the John mapping T . We obtain also 
two sets of points P n and Q n on the ellipses B and A, respectively. 

The mapping Ub has an inverse and generates a discrete dynamical system or, in 
other words, an action of the group ZonBas above for the John mapping. An orbit 
of this action is the set of the points P/., k = .... —1, 0, 1, 2, ... and Pfc = Ug^Pi. 
Now note that in a general case of a disposition the ellipses can be intersecting, then 
we must start from a point Q± on the ellipse A which is inside B and can determine 
the mappings Ub , Ua in the same way. In this case we can encounter with that 
the tangent straight line intersects the ellipse B only at one point, then we must 
consider this point as double. The ellipses can be tangent in one or two points, 
they can be tangent and intersecting simultaneously, they can be nonintersecting 
and lie one outside other, finally, they can be arbitrary irreducible conies. The 
way is the same. Note that one can consider the case the conic B is reducible, 
e.g., it is two nontangential different straight lines. Note more, that each projective 
transformation of the plane transforms a Poncelet mapping of conies into the same 
mapping of their immages, therefore we could restrict ourself to a case when one of 
conies is the unit circle. 

The first interesting problem is to describe these crossing points explicitly. It 
was solved by Jacobi and Chasles who showed that the sequences P n and Q n can be 
parameterized by the elliptic functions. The second problem, the so-called Poncelet 
porism or the big Poncelet theorem (see |13p. consists in showing that if a particular 
trajectory of action on conies is closed (i.e., if P/v = Pq for some N > 2) then this 
property does not depend on the choice of the initial point Q\ on the conic A (or on 
the choice of Pq). A modern treatment of this problem from an algebro-geometric 
point of view can be found, e.g., in |28j . Our approach is another. If we introduce 
the standard rational parameterizations of our conies (see below (|4.3|) . (14. 4|) ) then 
the parameters x and y of the points Q\,P\ prove to be connected by a polynomial 
equation F(x,y) = 0. In a generic situation F(x,y) is a polynomial of the exact 
degree 2 with respect to each variable. Indeed, for a nondegenerate situation a 
tangent to the conic A at a point x should intersect B in two distinct points, and, 
moreover, from a point y on B there are two distinct tangents to A. Conversely, for 
any polynomial F{x, y) of the degree 2'm. x and y it is possible to find two conies A 
and B parameterized as in (|4. 3[) and (|4.4[) . These arguments are sufficient to build 
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trajectories of the John dynamical system. Note that similar considerations were 
exploited in [25j in order to show that the tangent and intersection points belong 
to an elliptic curve. The authors of [25] also introduced two involutions I\ and I 2 
which are equivalent to our ones in the John T-algorithm for the Euler curve. In our 
approach we derive the curve F{x 1 y) = explicitly and study it. Note connections 
with Gelfand's question [38] and elastic billiard [58], [4Tj . 

4.3. Passage to the John mapping on a biquadratic curve. Find an explicit 
expression for the Poncelet mapping Ub- We introduce the standard rational pa- 
rameterization of an arbitrary conic [13j that can be found even if the conic is 
reducible. Assume that the conic A is described by the coordinates £o>£ij£2 of a 
two-dimensional projective plane, that is 

2 

i,3=0 

and, as it is well-known, the conic is irreducible iff the matrix A is nondegenerate. 
Note that we will describe vectors (contravariant tensors) by lower indexes and 
covectors (covariant tensors) by upper indexes because we often use indexes for 
exponents of degree. 

Corresponding affine coordinates will be denoted by £, 77: [£ : £1 : £2] = [1 : £ : T)], 
of course for £0 7^ 0. Then it is possible to find polynomials Eq(x), Ei(x), 23a (a;) with 
deg(Ei(x)) < 2 such that 

& = or £=fM, , = fM. (4 .3) 

E (x) E {x) 

Quite analogously, the conic B can be parameterized as 

G {y) G a {y) 

where Gj(y) are some other polynomials of the degrees not exceeding 2. Thus the 
value of parameter x completely characterizes the point on the conic A, and the 
value of y completely characterizes the point on the conic B. 
Let us consider more general case now. 

Lemma 1. Let curves A and B be given parametrically by (|4.3| . (|4.4|) with some 
smooth functions Ei, Gt and Eq ^ 0, Go ^ on each interval. The point P — G(y) 
lies on the tangent line to A at Q = E(x) iff F(x, y) = where 



F{x,y) 



E (x) E x {x) E 2 (x) 
E' (x) E' x {x) E' 2 {x) 
G (y) Gxiy) G 2 (y) 



(4.5) 



Proof. The affine tangent line L to the curve A at a point Q with the parameter 
x has the direction vector r = (4*, 42). Therefore the affine point P satisfies the 
equality 

OP(y) - OQ{x) = kr 

with some fceM and an origin O. The last equality is equivalent to the complanarity 
of vectors E = (E (x), E x {x), E 2 (x)), E' = (E' Q {x), E[ (x), E' 2 (x)) and G = (G (y), 
Gi(y),G 2 (y)) in the bundle space R 3 \ {0} of the projective fiber bundle M 3 \ 
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{0} -> MP 2 , where the prime denotes the derivative with respect to x. Indeed, the 
collinearity of the vectors OP(y) — OQ(x) 



E (x) 



Gi(y) 
G (y) 



E 2 (x) 
E a (x) 



G 2 (V) 
Go(y) 



and t means 
1 



Ex 
E 

Ei Gi 

Eq Go t 
i? 1 E\E[. 

■J— 



E, 



E, 



o 



Ei 
E 

E 2 G2 

Eq Go t 
E 2 _ E 2 E 
E ~ET 



E'q/Eq 



Ei/Eq 

E[/E 



E 2 /E 
g 2 /gq 
E' 2 /E{) 



E (x) 
E' (x) 
Go(y) 



Ei{x) 
E[(x) 
Gi(y) 



E 2 {x) 
E 2 (x) 
G 2 (y) 



Go(y)E (x) 



□ 



Return to our case of conies. We see F(x, y) has the form 

F(x,y) = M {x)G (y) + M x {x)G x {y) + M 2 (x)G 2 (y) (4.6) 

with polynomials Mi{y) defined as 

Mi(x) = e m (E' k (x)E l {x)-E k (x)E' l (x)), i,k,l = 0,1,2, (4.7) 

where e%ki is the completely antisymmetric tensor. One can easily check that 
deg(Mi(x)) < 2, therefore the curve F(x,y) = is a biquadratic curve of the 
form (TO)) . 

Note that the equality (|4.6p can be written as F(x,y) = ( M(x), G(y) ) with 
the scalar product (•, •) or F = (Al(x),G(y)} with the pairing (•, •) (in this case 
G e K 3 ,M e (K 3 )* ) and F = (E(x),E'(x),G(y)) with the mixed product in R 3 . 
Further, we introduce vectors x — colon(\^x^x 2 \y — colon(l,y 7 y 2 ) and matrixes 
E, G by the rules Ei(x) = Y^j^Eijxi = (Ex) i7 G = Gy, then the decomposition 

(|4.6|) can be written in the form F(x,y) = (Mx,Gy), where M(x) = Mx. This 
implies that F = (x,M*Gy) = (G* Mx,y). Comparing with the decomposition in 
the forms (|3.ip and (|3.2p we obtain 

A = G*M, B = M*G = A*, (4.8) 

where the matrix A is obtained in the same way as the matrix E above, moreover 
A = (an-) with the matrix from (|1.4j) . 

In a case of irreducible conies the matrix A is nondegenerate. Indeed, if the 
matrix A is degenerate then either the matrix M or the matrix G is degenerate 
by virtue of (|4.8[) . The degeneracy of G means a linear dependence of polynomials 
Gi(y) that is the conic B will be a straight line in this case. The degeneracy of M 
means that there exist constants cq 1 c\ 1 c 2 such that 

E Q (x) E 1 (x) E 2 (x) 
E' (x) E[{x) E' 2 {x) =0, 
c ci c 2 

i.e. a linear dependence of polynomials E^x) and the conic A will be a straight 
line. 

Return to the Poncelet construction. We have obtained the parameters x\ of 
the point Q\ and y\ of the point P\ satisfy the equation F{x\ 1 yi) = 0, with 
F from (14. 5[) . Note now that instead the point P\ we could write the point P 2 
with the parameter y 2 in the Poncelet construction and have the same equation 
F(xi, y 2 ) = 0. We obtain the first result: for any point given by a parameter value 
x\ on the conic A, the points with parameters yi and y 2 of the intersection points 
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of the tangent line L\ at x\ with B are determined as two roots of the quadratic 
equation: F(x\,y) — 0. Thus identifying a point with its parameter value we can 
say the Poncelet mapping Ub maps the point y\ into the point yi and hence it 
(more precisely, the mapping Is) coincides with the John mapping I\ from section 
2. It is analogously, the mapping U a maps the point x\ into the point Xi and 
hence it (the mapping I a) coincides with the John mapping I2, and the mapping 
U coincides with the John mapping T on the curve C (|1.4p . 
We have proved the following 

Proposition 20. Each pair of distinct irreducible real conies A, (|4.3[) and B, (|4.4p 
generates a biquadratic real polynomial F(x, y) of the form (|1.4p with nondegener- 
ate matrix A by means of (|4.5| such that the point y € B lies on the tangent at the 
point x € A if and only if F(x, y) — 0. The Poncelet mapping Ub gives us (and can 
be obtained from) the mapping U : C — > C that coincides with the John mapping 
T on the curve C (TO)) . 

Let us prove the following converse 

Proposition 21. For any biquadratic real polynomial F{x,y) of the form (|1.4p with 
nondegenerate matrix A and for every its decomposition (14. 6|) with given poly- 
nomials Gi(y), Mi(x) of the second order there exists an unique projective set of 
polynomials Ei(x) of the second order such that relations (|4.7|) (and (|4.5| ) hold, 
and, hence, we can relate with any such curve Fix, y) = and its decomposition 
(|4.6p a pair of conies A and B parameterized as in (|4.3p and (|4.4[) . 

Proof. Let us have a decomposition (|4.6p of a given biquadratic polynomial F . 
Then we have a parameterization of the conic B by means of projective coordinate 
Go,Gi,G2 and the first our aim is to find polynomials Eq,E\,E2 such that the 
equalities (|4.7p hold. In other words, we must find a polynomial solution of the 
system of ordinary differential equations 

Ex E' = M (4.9) 

with a known vector M = (Mq(x), M\{x), M2(x)), an unknown vector E = (Eq(x), 
Ei(x),E2(x)) and the vector product x in M 3 . One can interpret this equation as a 
problem to find a curve (more exactly, the tangent vector of a curve) if it is known 
its binormal. Now we will need the following 

Lemma 2. Consider the differential equation in K 3 

rxr' = b (4.10) 

with a known smooth vector function b — (bx(t), 62(^)1 ^(t)) and unknown vector 
function r = (ri(t), ^(i), r3(i)) depending on a real parameter t. If (b,b',b") < 
then the equation (|4.10p has not any smooth solution. If (b, b', b") > then a 
solution of the equation (|4.10p exists there and it is only r = ±(6, b', b")~ x / 2 b x 6'. 
If (&, 6', b") — on an interval then only four cases are possible with some parameter 
change s = s(t): 1) r = ro = const, 2) r = r^ + r±s, ro, 1*1 are constant vectors. 3) 
r = r$e s + r\e~ s and 4) r — ro coss + ri sins with the same r$,r\. 

Proof. Note first that the Jacobi determinant of the system (|4.10p of ordinary dif- 
ferential equations with respect to r 1 is identically equal to zero, so that a standard 
theory of differential equations systems does not work for this system. 
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I) . Let b x b' ^ on an interval. We observe that the solution must be only of 
the form r — v(t) b x b' because from (|4.10p we easily obtain: b ■ r = 0, b ■ r' = 0, 
so b' ■ r — 0. The scalar v is unknown still. Substituting such r in (|4.10p we obtain 
v 2 (b, b', b")b = b so that there exists no solution of (|4~T0|) when (b, 6', b") < 0. The 
equality (b, b' , 6") = on an interval gives b = 0, i.e. a contradiction. 

II) . Let now b x b' = on an interval. If b ^ and 6' ^ then b(t) = fi(t)bo 
with a scalar function fi, bo is a constant vector that means the curve r is a plane 
curve. Note that we examine simultaneously the case 6^0 and b 1 = (/i = 1). 
Substitute this 6 in (|4.10[) and choose another parameter t such that r x r' = bo- 
lt follows from this that r x r" = 0, so r" — v(i)r with a scalar function v. 
After a reparameterization t — > s we obtain r" — ±r, i.e. r — roe s + r\e~ s or 
r = ro cos s+r*i sin s with some constant vectors ro, ri of a plane which is orthogonal 
to bo. For ro ^ r\ these solutions are cases of a hyperbolic and an elliptic rotation 
of the plane pespectively, besides, there are also cases of a movement along a direct 
straight line (cases 6 = and ro = T\ ^ 0) and of a stationary state. □ 

Continuation of the proposition proof. We will apply the lemma to the case of 
polynomial vectors of the second order b = M and r = E. It is easily seen that for 
such M the equality M x M' = implies only trivial cases 1) or 2) of the lemma, so 
that M and v 2 = (M, M' , M"). The last mixed product doesn't depend 

on x. Indeed, 

-f-(M, M' , M") = (M\ M', M") + (M, M", M") + (M, M', M'") = 
ax 

as the vector M consists of quadratic polynomials. Note that we can change a sign 
of the mixed product (M, M', M") by a change of the sign of M. Thus, for given Mj 
from the decomposition (|4.6p there exists the scalar v such that v 2 = [M, M' , M"), 
hence the components of the derived vector E = vM x M' are polynomials of the 
second order (and the mixed product (E, E' , E") doesn't depend on x). □ 

4.4. Projective invariance of the biquadratic curve. One can consider the 
vector E x E 1 as a covector E* of a dual space (R 3 )* with components E^(x) = 
Y^i=o A t: *Ei(x), where the matrix A is from (|4.2[) . more exactly, E* is proportional 
to ExE'. Indeed, by definition we at once obtain (E*(x),E(x)) = and (E*,E') = 
0. The covector field E*(x) describes a conic in the dual projective space, the field 
of tangent lines to A (see, e.g., [2"T]): 

2 

.4,., /'. '/•■ ' (I (4.11) 

i,3=0 

with an inverse matrix (Aij) = (A 11 ) . This point of view gives us a representation 

F(x,y) = (E*(x),G(y)) (4.12) 

and help to understand why the equality E x E' = M with quadratic polynomial 
vector E implies M x M' = E in the projective sense as E** = E, but it doesn't 
too help for the case of a curve A of order larger than 2. 

Proposition 22. For each pair of different irreducible real conies A, (|4.3[) and B, 
(|4.4[) , for any projective transformation L of the projective plane the images LA 
and LB of conies A and B generate the same equation F{x, y) = and the curve 
(|1.4p . Two different pairs of conies with the same curve are connected by means of a 
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projective transformation A — > PlA, B — > Pl-B and give different decompositions 

Proof, fn order to derive the first result let remember that an arbitrary nondenerate 
projective transformation Pl of the affine plane can be obtained from a nondenerate 
linear transformation L of the bundle space R 3 \ {0} of the projective fiber bundle. 
Such linear transformation L gives a number factor det L for F in the formula (14.51) . 
Thus the arbitrary nondenerate projective transformation of the affine plane and 
conies A and B don't change the equation (|f .4[) . 

Further, given decomposition (|4.6[) can be written in the forms (13. f I) and (|3.2|) 
and generates the nondegenerate matrixes M, G and then A, B my means of (|4.8|) . 
Some another decomposition gives matrixes Mi, Gi and the same F and A: A = 
G*M = A a = GJATi, so that the linear transformations L = GiG -1 = (M^ 1 )* M* 
and MiM -1 = (G 1 " 1 )*G* translate the conies A, B into conies Ai = LA, B x = LB 
respectively as a projective transformations. Note from (|4. 2|) that under trans- 
formation L the matrix A passes into L~ 1 *AL~ 1 and the matrix M passes into 
L- l *M. □ 

Thus, the biquadratic curve F(x, y) = depends only on a projective class of 
the conies pair. But the conies can be in one of the following generic dispositions 
(the order is chusen for a correspondence to proposition 1 1 5|) : 

0) the conic B lies inside of A, strictly or not, i.e. no tangent line to A cuts B. 

1) the conic A lies strictly (i.e. without contact) inside of B so that each tangent 
line to A cuts B in two different points. 

2) the conic A cuts the conic B in four points and there exist a straight line that 
has no common point with A and B. 

3) the conic A intersects the conic B at only two points without contacts. 

4) the conic A lies strictly outside of B and the conic B lies outside of A. 

5) the conic A cuts the conic B in four points and each straight line has a common 
point with A or B (there is no common tangent line, a hyperbola and an ellipse). 

Now we would like to clarify how one could know by F what case occurs. In 
order to clarify it let's take the following observations. If the conic A intersects 
B in a point P then from the point P E B there is only a tangent line to A, so 
that in the Poncelet construction there is only a point (x, y) with the parameter 
y corresponding to P G B that is the point (x, y) is an y-vertex of G, i.e. an 
extreme point along direction of real axis y. Remember that by the proposition 
HH/or generic cases l)-5) the roots 0/1)2(2/1) and only they are y-coordinates of 
y-vertices. Each such vertex corresponds to a point of intersection A (~l B. 

In the same way a common tangent line to conies A and B in the plane (£, 77) 
gives us an rc-vertex (x2,y2) in the plane {x,y) and we have: if ^2(^2) = then 
2/2 = -A (x 2 )/Ax(x 2 ) otherwise D 1 (x 2 ) = A\(x 2 ) - ^Aq(x 2 )A 2 (x 2 ) = 0, y 2 = 
— Ai(x 2 )/2A 2 (x 2 ). And for generic cases l)-5) real roots of Di(x 2 ) and only they 
are x-coordinates of x-vertices. Each such vertex corresponds to a common tangent 
line. 

Let us remember that the statement |H] says the invariants g 2 and gs of polyno- 
mials Di and D 2 are the same. Note that the case A < (i.e. k 2 < 0) corresponds 
only to the above case 5) because the polynomial D\ has two real and two complex 
roots if and only if A < 0, therefore if and only if the polynomial D\ is the same. 
We obtain the following 
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Proposition 23. In the case of disposition 0) each of equations Di(xi) = and 
£>2(yi) = has no real solution and the curve is not real. In the case 1) each of 
equations Di(x±) — 0,1)2(2/1) = has no real solution, the curve C is real and 
it has no a vertex. In the case 2) each of equations Di(x\) = 0, £'2(2/1) = has 
four real solutions, the curve C has four x-vertices and four y-vertices. In the case 

3) each of equations Di(x\) = 0,1)2(3/1) = has two real solutions, the curve C 
has two a;- vertices and two y- vertices. In the case 4) we have four common points, 
so that the equation D\(xi) = has four real solutions, the equation 1)2(2/1) = 
has no real solution, the curve C has no y-vertex and four ^-vertices. In the case 
5) we have four common tangent lines, so that the equation D\(xi) = has no 
real solution, the equation 1)2(2/1) = has four real solutions, the curve C has no 
2;- vertex and four y- vertices. 

Any point of contact gives a point (xi,yi) which is both a;- vertex and y- vertex, 
i.e. Di(xi) = 1)2(2/1) = 0, the conditions (|3.18p of singular point are fulfilled, 
and according to the proposition [7] either the curve is irreducible or x\ and y\ are 
multiple zeros of the discriminants D\{x) and 1)2(2/) respectively and this singular 
point is unique. In both cases we observe either the case III of the John's list of 
dynamical system behaviors and a breakdown of smoothness of the curve or a case 
of degeneration. In this work we do not consider these cases. Some examinations 
are in the work [40] . 

Remark more that the statement [7] and propositionfTSIof the subsection l3.4l implv 

Proposition 24. For any two conies that are in one of generic dispositions l)-5) there 
are real linear- fractional changes R\ , R2 of real parameters x — R\ (x) , y — R2 (y) 
such that the corresponding reduced biquadratic curve C will be having one of 
the canonical form (i), (ii) or (iii) of the subsection 13.41 Here the disposition 1) 
corresponds to the case (iii), —2 < b < 2 ; the disposition 2) corresponds to the case 
(iii), b > 2; the disposition 3) corresponds to the case (iii), b < — 2; the disposition 

4) corresponds to the case (i), c = 1, b > 2; the disposition 5) corresponds to the 
form (ii), c = — 1; the disposition 0) corresponds to the case (i), a > 0, c = 1 with 
the condition of disappearence b < 2 i.e. b 2 < (a + l) 2 . 

4.5. Periodicity of Poncelet mappings. Now we can apply the John mapping 
T in order to construct two sequences of points x„ and y n on conies A and B. 
According to proposition [20] of the subsection 14.31 the Poncelet mapping Ub gives 
us (and can be obtained from) the mapping U : C — -> C that coincides with the 
John mapping T on the curve C (|1.4[) . Therefore the periodic trajectories in the 
Poncelet problem correspond to closure orbits of the John dynamical system. Thus 
we obtain main result 

Proposition 25. The Poncelet problem in generic setting is periodic if and only if 
the John mapping on corresponding biquadratic curve (|4.5p is periodic and then 
their periods are coincide. 

Explicit answers for the John mapping are given in propositions [15l [T8l [HI We 
obtain 

Proposition 26. Dispositions l)-5) of previous subsection correspond to cases l)-5) 
of the proposition [HI The conditions (|3 .41 [1 . (|3.49p - (|3.5ip are criterions of period- 
icity of the Poncelet mapping U in corresponding cases. 
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Note that last proposition implies the statement of the big Poncelet theorem. 

Since the Poncelet problem is invariant under the arbitrary projective transfor- 
mation of the plane (£, 77) we can reduce the conies A and B to some simple shapes. 
Consider several possibilities. 

(i) If we reduce A and B to concentric quadrics determined by the equations 

flat + rf/bx = 1, C 2 /«2+?7 2 /&2 = 1. 
In this case we have the parameterization 

£ = 2a\ /2 y/{l + y 2 ), „ = b\ /2 {l y 2 )/{l + y 2 ) 

for the conic A and 

£ = 2a\' 2 x/{l + x 2 ), 7? = b\ /2 {l - x 2 )/(l + x 2 ) 

for the conic B. It is easily verified that the polynomial Z(x,y) defined by (14. 5|) is 
reduced in this case to the simplest Euler-Baxter form: 

Z = x 2 y 2 + 1 + a(x 2 + y 2 ) + 2bxy 

with complex coefficients. 

Hence, we have a simple form solution of the Poncelet problem in this case as 
above (see (|3.39[) ): 

x n = Vk sn(h(n + s );k), y n = Vk sn(h(n + s + 1/2); k). (4-13) 

Parameters a, and hi can still be specialized further. For instance, it is possible 
to choose 02 = hi — 1 reducing B to the unit circle. This choice corresponds to 
so-called Bertrand model of the Poncelet process [50]. If the second conic is an 
ellipse £_ 2 /a 2 + T] 2 /b 2 — 1, a\ > 1, b\ < 1 then the formula (14. 5|) gives us a bounded 
curve of the form 

x2y2 + l + h {x2 + y2) 4ai + 1 = Q 
1 — 0\ 1 — 0\ 

In this case the points x n are isomorphic to the godograph distribution of spins in 
the classical Xy-chain (see section FTT1 and [26]). Another possible choice a 2 — a\ = 
62 — b\ corresponds to the confocal quadrics. In this case the Poncelet problem is 
equivalent to the elastic billiard [55], [4"T] . 

(ii) Let us consider a possibility to reduce the conies B and A to two parabolas 
in the euclidean plane (£,f?). One of them, say B, can be fixed by the choice 
£ = x,r] — x 2 , whereas another parabola remains arbitrary: £ = Fi(y)/ F (y), rj = 
-p2(y)/-Po(j/)i where F (y) = (ay + b) 2 is a square of the linear function (this is a 
characteristic property of any parabola). Then we get 

Z(x,y) = x 2 F (y)~ 2xF 1 (y) + F 2 (y). (4.15) 

Performing additional projective (complex) transformation of the variable y we can 
fix polynomials -F\(y), F 2 (y), and F (y) in such a way that the polynomial Z(x,y) 
becomes symmetric in x, y. Then we can reduce Z(x, y) to the form 

Z = (xy + (x + y)y + g 2 /A) 2 - (x + y + y )(4xyy - 53), (4.16) 

where g 2l 53 are two remaining (arbitrary) independent parameters of the polyno- 
mial Z. 

In this case we obtain the Poncelet points parameterized by the Weierstrass 
function 

x n = Aip(h(n + s )) + A , y n = A\p(h(n + s + 1/2)) + A . 
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As we will see in Sect. 7, the biquadratic curve of similar form appears for the 
phase portrait of the elliptic solution of the Toda chain. 

Finally, we note that when Fq — const and F±(y) is a linear function in y then 
the parabola A has its axis parallel to that of the parabola B. But then the 
curve Z(x, y) = describes arbitrary conies in the coordinates x and y with absent 
x 2 y 2 ,x 2 y, and xy 2 terms. 

4.6. Cayley determinant criterion. Let A, B be arbitrary conies as described 
in previous section. Recall, that all tangents are passing through points of the 
conic A, whereas all vertices lie on conic B. (It is possible to assume that A is 
located inside the conic B). Let Ma and Mb are 3x3 matrices describing conies 
(i.e. corresponding quadratic forms) in projective coordinates xq, x±, xi- Le. if the 
conic A is the unit circle x 2 + y 2 = 1 and conic B is concentirc circle with the radius 
R then quadratic forms for A, B are 

(4.17) 



x 2 



x Q and 



2^2 



R'x, 



Corresponding matrices Ma, Mb are diagonal: 

M A = diag(l, 1, -1), M D = diag(l, 1, 
Compute the characteristic determinant 

F(z) = det(A-zB). 



-R 2 ). 



(4.18) 



Clearly, F(z) is a cubic polynomial. This polynomial is invariant with respect to 
any similarity transformation A —>■ S~ 1 MaS, B — ► S~ 1 MbS with a nondegenerate 
matrix S. As well known, for a pair of quadratic forms, in general, apart from 
degenerate cases, there exist a transformation reducing both forms simultaneously 
to a diagonal form. The roots Zi,i — 1,2,3 of the polynomial F(z) have simple 
meaning. If the matrix Mb is reduced to identity matrix (i.e. Mb = diag(l, 1, 1)), 
then Zi,Z2, Z3 are diagonal elements (eigenvalues) of the matrix A. 

The first step in the Cayley criterion is Taylor expansion of the square root of 
the polynomial F{z): 

hc„/ + ... (4.19) 



\J F(z) = c + c\Z + c 2 z + • • • + 
Then we compute the Hankel-type determinants: 



and 





C3 


c 4 


Cp+l 


p 


C4 


c 5 


Cp+2 




Cp+l 


Cp+2 • • 


C2p_l 






C3 


■ Cp+l 


= 

-L-Lp 


C3 


c 4 


Cp +2 




C P+1 


Cp+2 • 


C 2 p 



Then the Cayley criterion [13] . [28] is: 



, p = 2,3,4, 



(4.20) 



P= 1,2,3, 



(4.21) 



Statement 8. The trajectory of the Poncelet problem is periodic with the period N 



if and only if = 0, if N = 2p and 
of the Cayley criterion see [28] . 



0, if N = 2p + 1. For modern proof 
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Illustration. Take again the simplest case of two circles (|4.17[) with radii 1 and 
R. Then obviously, 

F(z) = (z - l) 2 (zR 2 - 1). (4.22) 

The first nontrivial Taylor coefficients of y/F(z) are C2 = R 2 (R 2 — 4)/8 and C3 = 
R 4 (R 2 — 2)/16. The case C2 = means R = 2 which corresponds to perfect triangle 
trajectory. The case C3 = means R = 2 1 / 2 which corresponds to square trajectory. 

5. The Pell-Abel equation 

5.1. Historical notes on the Pell- Abel equation and some equivalent prob- 
lems of analysis. The Pell equation 

P 2 - RQ 2 = L (5.1) 

is a well-known Diophantine equation where for given integer number R which isn't 
a square, one must find integer P,Q,L satisfying the equation. It is established 
that L.Euler linked it to Pell by misunderstanding. Actually the equation (|5.1[) 
arose and was being studied in works of Indian mathematician Brahmagupta as far 
back as 1000 years before Euler (628 year) and also in works of Euler's predecessors 
(among them is P.Fermat). Theory of the Pell equation is well-known ([15]), note 
only that the standard theory of the equation (|5.ip establishes a connection with 
continued fraction of the number \/R~. 

The equation (|5.ip in a ring of polynomials R[t] or C[t] of one variable with a 
constant L is called the Pell- Abel equation (it may occur the name " Abel equation" 
and also "Pell equation for polynomials"). This equations arose in Abel's work 
of 1826 year where he studied representation of the primitive J p(t)/ \J R(t) dt in 
elementary functions, here p, R are polynomials. N.H. Abel proved that if this 
primitive can be represented by logarithm and rational functions of t and ^/~R then 
one can find polynomials P, Q and a number A such that 



B= d t = A^ P + ^ Q . (5.2) 

Ir p - \[rq 

Here the degree of R is even: degi? = 2m, degp = m — 1 and p/A = 2P' /Q. A 
main point of Abel's considerations is that the polynomials P and Q satisfy the 
equation (|5.1[) with L — 1. And conversely, if the polynomials P and Q satisfy the 
equation (|5.ip with L = 1 then the equality (|5.2[) holds with p — 2P' /Q and A = 1. 

Thus, the solvability of the Pell- Abel equation plays a role a integrability criterion 
of Abel differential. It is well-known that later J. Liouville, V.V. Golubev and 
successors shown if the integral in the left-side part of (|5.2j) with some p and R can 
be written by elementary functions then the right-side part must have the view (15. 2p 
with some P and Q and the same p. Note also, Abel gave one more criterion for 
the representation (|5.2p . namely the formula (|5.2p holds iff polynomial continued 
fraction of the function \f~R is periodic. Thus, the solvability of the Pell-Abel 
equation plays also a role of a criterion of periodicity of continued fraction. 

Consider one more classical problem - Chcbyshev's problem on searching of a 

polynomial of least deviation. Let us have a system of I closed intervals of real axe 

1-1 

/ = [—1, 1] \ 1J (a,j, bj) and the following problem. One should find a polynomial 

of given degree n with unit leading coefficient that has a least deviation on the set 
/, i.e. find a minimum of the functional \\t n — P n -i(t)\\ C (jy A general polynomial 
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P n -i(t) runs a finite-dimensional lineal space and the problem can be interpreted 
as a problem to find an element in a finitc-dimcnsional lineal subspace of a Banach 
space which is nearest to a given. But because the space C{I) isn't reflexive such 
nearest element isn't obliged to exist. 

By a method going back to P.L. Chebyshev, one succeeded in proving the ex- 
istence of such minimal polynomial 2J. Further, if the polynomial P„_i is min- 
imal on the set I then, possible, it will also be minimal on some greater closed 

m — 1 

subset E C [—1,1]. Such most wide subset E = [—1,1] \ (J (aj,(3j) is called 

j=i 

n-right extension of the set / |51j . If now one take a polynomial R in the form 

m — 1 

R = (t 2 — 1) Y[ (t ~~ a j)(t ~ Pj) then it is proved to be that the solvability of the 
j=i 

Pell-Abel equation (|5.ip with unknowns P(t), Q(t), L — const > is equivalent to 
that the set E is n-right extension of the set I. In addition the polynomial P gives 
a solution of extremal problem, and the number \fL is a least deviation (minimum 
of the norm) [5T] . 

It is interesting that in this case the subset E is continuous spectrum (which 
is absolutely continuous and two-valued) of some bi-infinite Jacobi (3-diagonal) 
selfajoint real periodic matrix in the space I 2 if and only if the set E is n-right 
or that is equivalently if the Pell- Abel equation (|5.ip is solvable (see references in 
|51j). Remark that the Chebyshev's problem will be got for the case of one interval 
E = [—1, 1], R = (1 — t 2 ), and the Akhiezer's problem - for the case of two intervals 
(degi? = 4) and polynomial of the fourth order E = [— 1, a] U [b, 1], a < b, R — 
(1 — t 2 )(t — a)(t — b). Note that there are several solvability criterions for the Pell- 
Abel equation with a polynomial corresponding to the Akhiezer's problem, among 
them is well-known Zolotarev's porcupine (see, for instance, [13] )■ Note also recent 
Khrushchev's work [3H] on links to continued fractions. 

5.2. Connection between the Poncelet problem and the Pell-Abel equa- 
tion. The solvability of the Pell- Abel equation 

P{tf - R(t)Q(t) 2 = L 2 (5.3) 

with the polynomial R of even order, as we have noted above, has several equivalent 
formulations ([51]). In the work of V.A. Malyshev [43] there is a new solvability 
criterion given in an algebraic form that we will formulate for interesting for us case 
of the fourth order R = t 4 + dit 3 + ... + d±. Let us expand the root y/R in Laurent 
series in a neighborhood of infinity: 

oo 

vr= ( 5 - 4 ) 

j=— m 

and make up a determinant of Hankel type: 







c 2 ■ 


Ck 




r fe = 


c 2 


c 3 . 


Ck+i 


^ k — 1,2,3 




Ck 


Ck+i 


C 2 k-i 





Malyshev criterion ([13]) states: 
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Statement 9. The Pell- Abel equation (|5.ip with polynomial R of the fourth order 
has some polynomials P and Q of orders k + 2 and k as a solution if and only if 

r fe = o. 

Our observation is following. Let consider the Pell-Abel equation (|5-3[) with R 
of the fourth order and let Ai be one of roots of polynomial R(t). Make the shift 
of parameter t — > t + Ai (i.e. t = t + Ai). New Pell- Abel equation (|5.3|) will 
also be solvable. Now apply the Malyshev's criterion to new equation and let the 
coefficients Cj be coefficients of Laurent series of the root */ ' R{t + Ai). Then the 
equality = is necessary and sufficient for solvability of the equation (|5.3|) with 
given orders of polynomials P, Q. On the other hand, let in initial setting make 
change of variable s = l/(t — Ai), £ = Ai + 1/s. Then we obtain 

R(t) = ^(s)/s 4 (5.6) 

with a polynomial -F(s) of the third order and therefore y/~R = ^F(s)/s 2 . Note 
that the polynomial R can be restored from F by inverse transformation: 

R = t 4 F{l/t) = t( Zl t - \){z 2 t ~ l)(z 3 t - 1), (5.7) 

where Zi, Z2, z 3 — roots of the polynomial F of the third order. The polynomial F 
can be generated by formula F = det{A — zB) and matrixes A — diag(zi, z 2 , £3), 
B = diag(l, 1,-1) that are built by quadratic forms Z\x\ + z 2 x\ + z 3 Xq and x\ + 
x\ — Xq of conies 

z lX 2 + z 2 y 2 + z 3 = 0, x 2 +y 2 = l. (5.8) 

The expansion (|4. 19|) of the root \J F(s) gives us the expansion y/ F(s)/s 2 = 
cq/s 2 + ci/s + c 2 + ■ ■ ■ + c n s n + . . • which after inverse change can be written as 

y / R(t + Ai) = c t 2 + Cl t + c 2 + c 3 /t + h c n /t n ~ 2 + .... 

We see that Ci = c 3 ,C 2 — C4, ... and that the determinant (|4.20p coincides with the 
determinant (|5.5p . i.e. the Cayley criterion for the case of the even period with the 
polynomial F coincides with the Malyshev's criterion for the equation (|5.3p with 
the polynomial R. We have proved 

Proposition 27. Let the Poncelet problem for given conies A and B generate a 
polynomial F of the third order by formula (|4. 1 8[) and a polynomial R by formula 
(15. 6[) . Then the equation (|5.3[) with polynomial R of the fourth order is solvable 
if the Poncelet problem is periodic with an even period. Vice versa, for given real 
polynomial R of the fourth order without constant term one can build a pair of 
conies A and B so that the Poncelet problem (generated by polynomial F by means 
of formulae (|4. 18[) and (|5.6[) ) is periodic with an even period if the equation (|5.3[) 
with polynomial R is solvable. 

Now the proposition [TO] gives us solvability criterions for the Pell- Abel equation 
6. Criterions of uniqueness breakdown for the Dirichlet problem 

AND RlTT'S PROBLEM 

In this section we give criterions of solution uniqueness for the Dirichlet problem 
(|l.l[) , (|1.2p and show how to construct a system of some special solutions in case of 
non-unique the Dirichlet problem for the string equation. We describe a method 
that can be applied to a slightly more general class of curves than biquadratics. 
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We'll consider the Dirichlet problem (|l.lj) . (11.2j) in the classical setting (|2.2j) if 
the curve C is bounded and in modified setting (12. 3|) if the curve C is unbounded 
because we'll need an application of the proposition [TJ 

Remind that in subsection 13.51 we consider an even periodic meromorphic func- 
tion 4>{z) with some (complex) period T and a complex curve C C C 2 described 
parametrically as 

x{t) = <j>{t), y(t)=(p{t + s), teC (6.1) 
with the property (|3.35p . £ is a nonzero complex parameter. We considered also a 
real curve C = C DM 2 given by means of a contraction of the functions x, y on a 
set S <E C. We saw that the curve C is symmetric with respect to the line y = x, 
the complex John mapping T is equivalent to a shift of the parameter t by the 
step —2e and the periodicity condition of the complex John mapping has the form: 
2ne = mT with some integer n, m. 

We assume also (this is a very strong restriction, as we will soon see) that the 
function (f>(z) possesses a nontrivial multiplication property: 

<f>(nz) = R n {<f>{z)), n = 1,2,3,... (6.2) 

where R n (z) is a family of rational functions of the argument z (by definition, 
R\(z) = z, but for other n = 2,3, . . . the expression and even degree of R n {z) can 
be non-obvious). We will consider the condition 

Vx e K, R n (x) ^ oo (6.3) 

Consider the complex setting (|2.7j) of homogeneous Dirichlet problem. Now the 
sufficient condition of the uniqueness of the Dirichlet problem from proposition [2] 
requires the transitivity of T that implies that for each integers n and m we have 
2ne 7^ mT. 

Assume now that this condition isn't fulfilled, i.e. for some integers N, M we 
have the condition 

2Ne = MT (6.4) 
We show that in this case the problem is indeed non-unique and we construct 
explicitly a system of explicit solutions 

$„(x, y) = f n {x) + g n (y), n = 1, 2, . . . (6.5) 

for the string equation in the domain bounded by the curve C . 

As the rational functions R n (z) are assumed to be non zero, bounded on R 2 and 
non-constant then the function $(x,y) = f n [x) + g n (y) — R2un{x) — i?2njv(2/) is 
obviously nonzero and smooth in R 2 . So we should verify the Dirichlet boundary 
condition <&(x,y) = in all points of the curve C. Indeed, for any point t in the 
curve we have 

fn(x(t)) + g n {y{t)) = R 2 nN(x(t)) - R 2 nN(y(t)) - (f>(2nNt) - <f>(2nNt - 2nNe) = 0, 

where we exploited properties (|3.34[) . (|6.2p and (16. 4[) . Hence the function $ n (a;, y) 
is identically zero in all points on the curve C. 
We proved the following 

Proposition 28. Under the condition (|6.4p for the curve C we can choose 

f n (x) = i?2njv(x), n (j/) = - i?2nJV (y) , (6.6) 

where R n (z) are rational functions defined by (|6.2|) with property (|6.3j) . This will 
provide a non-zero solution (|6.5p of the Dirichlet problem in the complex setting 
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()2.7j) . If, in addition, the condition (|6.3|) is fulfilled then obtained solution satisfies 
the modificed setting (|2.3p and if furthermore the curve is bounded then we have a 
classical solution in interior of the curve C, i.e. in the setting (|2.2|) . 

We first illustrate how this theorem works in simplest case when the curve C is 
ellipse. In this case we can choose parameterization 

x(t) = cos(t), y(t) = cos(t + e) (6.7) 

with some parameter e. Note that the ellipse described by (|6.7p has semiaxis 
a = \[2 cos(e/2), b = v / 2sin(e/2) one of which is inclined by the angle n/4 with 
respect to axis x. 

Now the John criterion (|6.4p of uniqueness of solution is 

e ^ Nn/M. (6.8) 

Otherwise, assume that e = Nn/M for some positive integers N,M. The function 
cos z satisfies all needed conditions, moreover, we have 

cos(zn) = T„(cos(z)), n = 1, 2, . . . (6.9) 

where T n (x) are the Chebyshev polynomials of the first kind, hence the proposition 
[2H1 gives us explicit solutions in this case 

f n (x) = T 2Nn {x), g n {v) = -T 2Nn {y), n = 1,2,... (6.10) 

Thus, in case of nonuniqucness of the Dirichlct problem for ellipse we get a 
dcnumerable family of polynomial nonzero solutions inside the ellipse. 

Let C be the biquadratic curve. In beginning, let's consider the case (|3.42j) of 
subsection 13.51 a > 0, c = — 1. Here the curve is parameterized as 

x = <j>{t) = cn(t; k), y = 4>{t + ri) = cn(t + r/; k). 

An easy analysis shows the last curve will be real for real t and real rj. For the 
elliptic cosine there are well-known multiplication formulae 

cn(2mz) = T^(cn 2 z), cn((2m + l)z) = an z T^cn 2 z) 

with some rational functions T^(z),T 2 (z). Poles of the function cn(z) lie on the 
line Imt — K' therefore the condition (|6.3| holds, the curve C is bounded in this 
case and the proposition [28] states the condition (|3.49| is sufficient for uniqueness 
breakdown for the Dirichlet problem in classical setting. Necessity of the condition 
(|3.49[) follows from the proposition [1] (for details see [51)]). 

In similar ways all another cases can be considered (it done in works [l8].[llj.[56j). 
We formulate final result: 

Proposition 29. The conditions (|3.49[) - (|3.51[) are criterions of uniqueness breakdown 
for the Dirichlet problem (|l.ip . (|1.2p in corresponding cases (|3.42[) - (|3.48[) of the 
biquadratic curve in a canonical form. The solutions are understood in modificed 
setting (|2.3p of the Dirichlet problem and for the cases of bounded curve C — in 
classical setting (|2.2|) . Under the condition there is a denumerable set of linear 
independed analytical solutions of the homogeneous Dirichlet problem. 

Now we can consider the problem: what is the most general class of the even 
periodic functions (f>(z) with the property (|6.2p ? This problem has solved by Ritt 
in 1922 [49]. The list of possible such functions includes two generic cases: 

W = cco S S (z)+d with arbitrary a, b, c, d; 
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(ii) 4>(z) is generic even elliptic function of second degree, i.e. cf>(z) — ^^jjg , 
or, equivalently <fi(z) — n ^Z^aft+e^) w ith arbitrary parameters a, b, c, d, k, e%, e^- 
Parameters 32,33 of the Weierstrass functions are arbitrary (equivalently, periods 
2o;i,2a;2 are arbitrary). 

There are also 3 exceptional cases connected with elliptic functions with certain 
restrictions upon their parameters: 

(iii) 4>{z) = (z)+d w ith arbitrary a, b, c, d but with restriction 33 = and 32 
arbitrary real. In this case the fundamental parallelogram of periods 2lo2,2lo3 is a 
simple square. Such case corresponds to so-called lemniscatic elliptic functions. 

(iv) <f>{z) = (z]+d w ith arbitrary a, 6, c, d but with restriction <?2 = and 33 
arbitrary real. In this case the period lattice has the hexagonal symmetry. Such 
case corresponds to so-called equianharmonic elliptic functions. 

For the equianharmonic case Ritt found one more case when 4>(z) is a linear 
rational combination of the Weierstrass function p'(z), but in this case the function 
4>{z) is not even and we can omit this case. 

The case (ii) corresponds to our problem in biquadratic. Thus, in the case of 
nonuniqueness we have the solution 

f n (x) = B.2nN{x), g n (u) = ~R2nN (y), (6-11) 

where R n (z) are rational functions defined as (16. 2|) . Note that in contrast to the 
trigonometric case (i.e. <f>(z) — cos(z)) there are no simple explicit expressions for 
R n (z) for arbitrary n. 

What about cases (iii) and (iv)? In these cases we have that the curve C is 
described by an algebraic equation of degree greater than 4. Consider, e.g. the case 
(iii) and assume for simplicity that <j>{z) — p 2 (z). We are seeking an algebraic equa- 
tion f(x,y) — for the curve C described parametrically as x(t) = p 2 (<), y(t) = 
p 2 (t + e). This equation can be easily found from addition theorem for the Weier- 
strass function. We will not write down it explicitly, let us mention only that degree 
of this equation is 8. Analogously, for the case (iv) we obtain the algebraic curve 
describing by a polynomial F(x,y) of a total degree 12. But we should note that 
the John condition of intersection at no more than two points with vertical and 
horizontal lines is not satisfying and here we do not have a sufficient condition of 
solution uniqueness for boundary value problems. 



7. Related problems of mathematical physics 

7.1. Classical Heisenberg XY spin chain. There is an interesting relation with 
the classical Heisenberg XY spin chain [26] which is a system of 2-dimensional unit 
vectors ("spins") r n = (x n ;y n ), \f n \ = 1 with energy of interaction 

N-l 

E= ^(r„, Jf n+1 ), (7.1) 

n=0 

where J = diag(Ji, J2) is a diagonal 2x2 matrix. The problem is to find static solu- 
tions providing local extremum to the energy E. As was shown in [26) this problem 
is equivalent to finding solutions of the systems of non-linear vector equations: 



(f n , J(r n -i + f„+i)) =0, n = 1,2, . . . ,N - 1. 



(7.2) 
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Below we would like to select cases of closed chain when r*o = r/v for some N. In 
what follows we will assume that J± = 1 , J2 = j > 1 . Among all solutions of (|7.2j) 
we choose so-called regular ones [55] with the condition 

r n -i + f n+ i ^ 0. (7.3) 

Then it is possible to show that the scalar product 

W=(r n ,Jr n+1 ) (7.4) 

doesn't depend on n and hence can be considered as an integral of the system 
(|7.2p . It is enough to construct general regular solutions [26] which are of two 
types. Choice of solution depend on the value of integral W. If \W\ < l/j then 

x n = cn(q(n - 9);k), y n = sn(q(n - 6);k), (7.5) 

where parameters k, q can be found from 

dn(q;k) = l/j, fc 2 = I—^J. (7.6) 

If l/j < \W\ < 1 then 

x„ = dn(q(n - 9);k), y n — k sn(q(n - 9); k), (7.7) 



where 



U2 



cn(q;k) = l/j, k z = — . (7.8) 

l-J 

In both cases parameter 9 is arbitrary real number depending on initial condition. 
If the chain is periodic then we have 

qN = AKmi + 2iK'm 2 (7.9) 

which coincide with (|3.40p . The reason for such coincidence is the following. 
Consider the equation of the integral 

x n x„+i +j~ 1 y n y n +i = W (7.10) 

with fixed value W. The equation (|7.10|) can be reduced to the algebraic form (|3.28|) 
by standard substitution (stereographic projection from the unit circle to the line): 

_ 1 - m„ _ 2u n 

1 + < 1 + < 

It is easily seen that the variables u„,u„+i lie on the Euler-Baxter biquadratic 
curve 

u l u n+i + 1 + a{u 2 n + u 2 n+1 ) + bu n u n+ i = 
with parameters a, b simply related to the "physical" parameters j, W. 

Then it is easily verified that finding solutions (step-by-step) of equations (|7.2[) 
for the regular solutions are equivalent to finding points M2, M3, . . . Mm for the 
John algorithm. Note that arguments similar to ones in the formula (|3ld|) give us 
that the parameters q and 9 must be real so that the condition (|7.9j) takes the form 

-5L = (7.11) 

Thus static regular solutions for the closed finite classical Heisenberg XF-chain 
are equivalent to periodic solutions of the John algorithm for the Euler-Baxter 
biquadratic curve, or, equivalcntly, to finding periodicity condition of the Poncelet 
process for the unit circle x 2 + y 2 = 1 and a concentric conic x 2 /£i + y 2 /^2 = 1- 
This choice of conies corresponds to the Bertrand model of the Poncelet process 
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[50j . This means that there is an equivalence between static periodic solutions of 
the Heisenberg XF-chain and the Bertrand model of the Poncelet process. 

7.2. Elliptic solutions of the Toda chain and biquadratic curve. Now we 

consider the Toda chain that is a discrete dynamical system consisting of two sets 
u n {t), b n (t) of complex variables depending on continuous parameter t and discrete 
parameter n — 0, ±1, ±2, .... The equations of the motion are [57] 

b n = u n+1 -u n , u n = u n (b n - b n _x) (7.12) 

The Toda chain is one of the most simple and famous completely integrable discrete 
dynamical systems (for history of this model and review of different approaches see 
[57|). Among numerous explicit solutions there are so-called "elliptic waves" con- 
structed firstly by Toda himself [57]- We give the Toda elliptic solutions in some- 
what different form which is more convenient for applications (Toda used Jacobi 
elliptic functions whereas we exploit the Weierstrass functions). 

Proposition 30. Elliptic solution for the unrestricted Toda chain can be presented 
in the from 

b n = w(c(ujt — p(n + 1) + r) — ((ut — pn + r)J — A, 

U n = ^ 2 (p(p) - pM - P n + r )j = 



,2 



a(ivt — p(n + 1) + r) o~{ujt — p(n — 1) + r) 



(7.13) 



a 2 (p) a 2 (tot — pn + r) 

with arbitrary parameters u>,p, r, A and arbitrary invariants <?2, 33- Here p(z), cr(z), 
((z) are standard Weierstrass functions defined as in [50] , 

Proof. The Toda chain equations (|7.12p are verified directly by substitution using 
well known formulas [60j : 

d >/ s , s d , , . . , . , . a(y — x)a(x + y) 

— C(« = -p( z ) ; ~r !g = CW, p(z) - P(y) = yy 2I \. u ' 

dz az <7 z {x)<j z (y) 

Strictly speaking, the variable b n (t) is inessential - it can be eliminated from the 
system (|T. 12[) . Thus only variable u n (t) is considered as a "true" Toda chain variable 
(for details see [ST]). 

Now we can construct the phase portrait for this variable. By the phase portrait 
we will assume the plot constructed on the plane x, y if one take points Po, P±, P2, . ■ ■ 
with the coordinates Pq — (uo,ui), Pi — (ui, 112), . ■ ■ , P n = {u n ,u n+ i). The 
phase portrait is an indicator of integrability: if the system is integrable the points 
Pi,i = 0,1,... fill some smooth curve. Otherwise these points are distributed 
quasi-stochastically (so-called "stochastic web" [3T]). In our case the variable u n (t) 
is given explicitly by (|7.13[) . u n +\(t) = u n (t — p/u>) and hence the points P n — 
(u n ,u n +i) fill the symmetric biquadratic curve in its canonical form (|3.24[) by the 
proposition[9] The complex John algorithm in this biquadratic curve is equivalent to 
passing from the point P n to point P n +i and then to point P n +2- Note that the time 
parameter t describes a smooth motion along this curve (Hamiltonian dynamics), 
whereas the John algorithm describes a "discrete motion". The periodic case is 
defined by a positive integer N such that u n (t) = u n+ N(t). From (|7.13|) it is seen 
that periodicity condition is equivalent to 

pN = 2m\U)x + 2m 2 w 2 (7.14) 
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with some integers mi,m2. An important property of the elliptic Toda solutions 
is that the periodicity property (|7.14|1 doesn't depend on the time parameter t, i.e. 
if periodic condition takes place for one value of t then it also holds for all others 
values of t. In terms of the (complex) John algorithm this means that a period of 
a point following the John mapping doesn't depend on choice of this initial point 
Pq = i u o(t) t ui(t)) on the curve. 

Comparing obtained formulas for u n (t) with (|4.15p and (|4. 16[1 we can say that the 
elliptic solutions of the Toda chain give an interesting illustration of the Poncelet 
theorem: periodic solutions (|7.14[) of the Toda chain correspond to periodicity of 
the Poncelet process on two parabolas. 



7.3. Elliptic grids in the theory of the rational interpolation. In this section 
we describe briefly an interesting connection of the Poncelet problem with theory 
od so-called admissible grids for biorthogonal rational functions. This subject is 
important in theory of rational interpolations. For further details and relations 
with theory of special functions see, e.g. [55], [53], [H], [53J, [54] , 

Let us consider the set of rational functions R n (x) of the order [n/n], which 
means that R n {x) are given by ratios of two n-th degree polynomials in x. In what 
follows it is assumed that all rational functions R n {x) have only simple poles. 

We take ai, a-z, . . . , a n as n distinct prescribed positions of the poles of R n {x). 
Then R n (x) can be written as a sum of partial fractions 



(n) 



i?n(z)=4 n) +E^ ( ? - 15 ) 

%— l 

with the coefficients t\ , i = 1,2, ... ,n, playing the role of residues of R n {x) at 
the poles ai. The coefficient can be interpreted as lim^—^ R n (x). 

Let x(s) be a "grid" , i.e. a function in the argument s. We would like to construct 
so-called lowering operator "D x r s ) in the space of rational functions R n {x) defined 
on this grid in the following way. We take as a definition of the lowering operator 
T^x(s) a divided difference operator in the parameterizing variable s, which obeys 
the following properties: 

i) the grid x{s) is a meromorphic function of s £ C which is invertible in some 
domain of the complex plane; 

(ii) for any function F(s) one has 

V x{s) F(s) = X (s)(F(s + l)-F(s)), 

where x( s ) is some function to be determined; 

(iii) T> x ( s )Ri{x) = const, where Ri(x) is an arbitrary rational function of the 
order [1/1] with the only pole at x = ot\] 

(iv) the operator T> x ^ transforms any rational function R n {x) with the pre- 
scribed poles a\,...,a n to a rational function R n -\(y(s)) of the adjacent grid y(s) 
with some other poles (3\, . . . , f3 n -x; 

(v) the operator X>w s ) is "transitive" : for any nonnegative integer j the operator 

V [J ) > defined as 

x(s) 

V%F(s)= Xj (s)(F(s + l)-F(s)) 
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with some function Xj ( s ) transforms any rational function R n {x) with the prescribed 
poles Q!j+i, ■ ■ ■ , Oj+n to a rational function i?„_i(j/(s)) with one and the same 
adjacent grid y(s) and the sequence of poles /3j+i, . . . , (3j+ n -i, 

(vi) we assume that the poles are nondegenerate: there are infinitely many dis- 
tinct values of a n and j3 n and a n 7^ a n +i, a n+ 2 and, similarly, j3 n 7^ /3„+i, (3 n +2 for 
all n. 

An important restriction is the condition of independence of 'D x ( s ) on the order n 
of a rational function. The problem is to deduce the functions Xj( s ) an d x ( s )>y( s ) 
from the given set of requirements. 

From the properties (i)-(iii) we easily find 

xis)= ( \ ?_v l 

\x(s + 1) — oti x(s) — OL\ ) 

'X (7.16) 



x(s) — x(s + 1) 

The function x( s ) is defined up to an inessential constant multiplier. Note that 
from (ii) we have V x ^Rq(z) = 0. 

The most non-trivial problem in the construction of the operator T> x ^ s \ consists 
in establishing the properties (iv)-(v). 

This problem was solved in [SI]. It appears that the grid x(s) as well as the 
grids y(s),a s , f3 s should belong to the class of the elliptic grids. This means that 
they should satisfy the biquadratic equation 

A 1 x 2 {s + l)x 2 (s) + A 2 x 2 (s)x{s + 1) + A 3 x 2 (s + l)x(s) + A 4 x(s + l)x(s) + 

A 5 x 2 (s + 1) + A 6 x 2 (s) + A 7 x(s + 1) + A s x(s) + A 9 = (7.17) 

with some constants Aj,i = l,2,...,9. As we already know this equation can be 
parameterized in terms of the elliptic functions and we thus arrive at elliptic grids 
x(s). Thus the elliptic grids x(s), y(s) are the most general ones to provide existence 
of the lowering operator for rational functions. 

In the theory of the rational (sometimes called the Cauchy-Jacobi or Pade) 
interpolation these grids appear naturally for some class of self-similar solutions. 

The Cauchy-Jacobi interpolation problem (C JIP) for the sequence Yj of (nonzero) 
complex numbers can be formulated as follows [S], [44j . Given two nonnegative in- 
tegers n, m, choose a system of (distinct) points Xj,j = 0, 1, . . . , n + m on the 
complex plane. We are seeking polynomials Q m (x; n), P n (x; m) of degrees m and 
n correspondingly such that 

Q (x 3 ;n) n + m 

F n (Xj;m) 

(in our notation we stress, e.g. that polynomial Q m (x;n), being degree m in x, 
depends on n as a parameter). 

It can happens that solution of the CJIP doesn't exist. In this case it is reasonable 
to consider a modified CJIP: 

YjP n (xj-,m) - Q m {xj-,n) =0, j = 0, 1, . . . ,n + m, (7-19) 

where polynomials P n (x; m), Q m (x; n) can be now unrestricted. The problem (|7.19p 
always has a nontrivial solution. In exceptional case, if the system (|7.18[) has no 
solutions, some zeroes of polynomials P n {z;m) and Q m (z;n) coincide with inter- 
polated points x s - Such points, in this case, are called unattainable [44] . 
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The CJIP is called normal if polynomials Q m (x; n), P n (x; m) exist for all values 
of m,n = 0,1,... and polynomials Q m (x; n), P n (x; m) have no common zeroes. 
This means, in particulary, that polynomials Q ' m (x; n) , P n (x; m) have no roots, 
coinciding with interpolation points, i.e. 

Q m (xf, n) ^ 0, P n (xj; m) ^ 0, j = 0,1,. ..n + m (7.20) 

In a special case when there exists an analytic function f(z) of complex variable such 
that f{xj) = Yj the corresponding CJIP is called multipoint Pade approximation 
problem [8]. 

It is possible to show that the Cauchy-Jacobi interpolation problem is equivalent 
to theory of biorthogonal rational functions |64j . The elliptic solutions (obtained 
hrst in [52]) of this problem appear naturally if the interpolated function f(z) satisfy 
the so-called discrete Riccati equation [42]. Geometric interpretation of obtained 
elliptic grids and their connection with the Poncelet problem can be found in [32] 
and [S3]. 

Note that if terms of degree > 2 are absent in (|7.17[) (i.e. A\ = = A3 = 0) 
then corresponding grid is degenerated to so-called Askey- Wilson grid [7] which is 
the most general grid for orthogonal polynomials satisfying a linear second-order 
difference equation 61j. From geometrical point of view the Askey- Wilson grids 
correspond to the John algorithm for the second-degree curves (i.e. ellipsis, hyper- 
bola or parabola) [42] . 
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